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Abstract 



Possible deviations from a radiation-dominated evolution, occurring prior to the synthesis of 
light nuclei, impacted on the spectral energy density of high-frequency gravitons. For a systematic 
scrutiny of this situation, the ACDM paradigm must be complemented by (at least two) physical 
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parameters describing, respectively, a threshold frequency and a slope. The supplementary fre- 

1 

quency scale sets the lower border of a high-frequency domain where the spectral energy grows 
with a slope which depends, predominantly, upon the total sound speed of the plasma right after 
inflation. While the infra-red region of the graviton energy spectrum is nearly scale-invariant, 
the expected signals for typical frequencies larger than 0.01 nHz are hereby analyzed in a model- 
independent framework by requiring that the total sound speed of the post-inflationary plasma 
be smaller than the speed of light. Current (e.g. low- frequency) upper limits on the tensor power 
spectra (determined from the combined analysis of the three large-scale data sets) are shown to 
be compatible with a detectable signal in the frequency range of wide-band interferometers. In 
the present context, the scrutiny of the early evolution of the sound speed of the plasma can then 
be mapped onto a reliable strategy of parameter extraction including not only the well established 
cosmological observables but also the forthcoming data from wide band interferometers. 
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1 The general famework 

Cosmological observations rely on three pivotal data sets, i.e. the Cosmic Microwave Background 
(CMB) data, the determinations of the matter power spectrum from galaxy surveys and the supernova 
light curve observations. The large-scale measurements are inextricably bound to the model used to 
interpret the data. Consequently the three aforementioned data sets are jointly analyzed in terms 
of a standard scenario which is often dubbed ACDM paradigm, where A qualifies the dark energy 
component and CDM denotes the cold dark matter component. While all the current observations 
are based, directly or indirectly, on the electromagnetic spectrum, there is the hope, in the future, 
that the electromagnetic observations could be complemented by the analysis of the spectrum of the 
relic gravitons which have been produced both in the context of the ACDM paradigm as well as in 
other related contexts. The problem is, therefore, twofold: on the one hand reliable estimates of the 
spectrum of the relic gravitons arising in the ACDM paradigm are needed. On the other hand it 
will be important to analyze other complementary scenarios. The purpose of the present paper is to 
address both issues in quantitative terms. The relic graviton background produced in the context of 
the ACDM paradigm is expected to be rather minute and undetectable by wide-band interferometers 
[H 121 [3J H] in one of their future realizations. There are, however, extensions of the ACDM scenario 
where the signal potentially detectable by wide-band interferometers is much larger than in the current 
paradigm. 

The recent WMAP 5-yr data [HI El 0, El [9] set quite stringent bounds on the amplitude of the relic 
graviton spectral energy density for typical frequency scales q 

k 
u p = -£■ = 3.092 x KT 18 Hz = 3.092 aHz, (1.1) 

where, according to the prefixes of the international system of units, 1 aHz = 10 -18 Hz. The wavenum- 
ber fc p is also called sometimes pivot scalqj since it is customary, in the experimental analyses of CMB 
data [8, 9), to assign the amplitude of the scalar and tensor modes exactly at k p . In the same units 
of Eq. (II. ip the typical frequency interval potentially accessible to the observations of the wide-band 
interferometers ranges between few Hz and 10 kHz with a peak of sensitivity around 100Hz. 

The logic followed in the present investigation will be to compute as accurately as possible the relic 
graviton spectra in terms of the parameters of the putative ACDM paradigm. The CMB data will 
then be used to enforce the normalization of the spectra at u p . This will allow for the estimate of the 
spectral energy density at the frequency explored by wide-band interferometers (i.e., approximately, 
100 Hz) not only in the case of the ACDM paradigm but also in the context of its extensions. The 
latter extensions will be examined on the basis of their plausibility, e.g. the models which are already 
incompatible (or barely compatible) with CMB observations will not be analyzed and the attention 
will be focussed on those scenarios which are not ruled out by (current) large-scale observations and 
which may lead to a potentially large signal at the wide-band interferometer scale. In the present 



2 Natural units % = c = k-g, = 1 will be consistently adopted all along the present investigation. 
3 The pivot wavenumber k p corresponds to an effective multipole £ e g ~ 30. 



introductory section, after a general discussion of the typical frequencies of the graviton spectrum, 
the present status of CMB observations and wide-band interferometers observations will be swiftly 
discussed. Specific attention will be paid to those quantitative aspects which are germane to our 
theme, i.e. the stochastic backgrounds of relic gravitons. Some of the concepts introduced here will 
also be more specifically addressed in the forthcoming sections. At the end of this introduction, the 
purposes of the present investigation will be more specifically outlined. 

1.1 Typical frequencies of the problem 

To compare frequencies it is mandatory to specify the background and the appropriate conventions 
on the normalization of the scale factor. Consistently with the ACDM paradigm, the background 
geometry will be taken to be conformally flat, i.e. 

ds 2 = g v dx' J 'dx 1 ' = a 2 (r)[dT 2 — dx 2 ], g = a 2 ^)!]^, (1.2) 

where 77^ is the Minkowski metric with signature mostly minus, i.e. (+,—,—,—). The scale factor 
at the present time will be normalized to unity, i.e. ao = 1. Within the latter convention, the 
comoving frequencies (or wavelengths) coincide, at the present time, with the physical frequencies 
(or wavelengths). The (conformal) time derivative of the logarithm of the scale factor will be used 
throughout the script and it is defined as 

n = a - = d -^, (i.3) 

a dr 

Note that, in Eq. (jl.3p . the prime denotes a derivation with respect to the conformal time coordinate r: 
this notation will be consistently enforced in the whole investigation. The evolution of the background 
can be expressed in terms of 7i and Ji! and it is given by: 

3H 2 = a 2 4p t , (1.4) 

2(H 2 -H')=a 2 eUp t + Pt ), (1.5) 

p' t + 3H(p t + Pt ) = 0, (1.6) 

where p t and p t denote, respectively, the total energy density and the total pressure of the plasma. 

The frequency of Eq. (11.11) can be usefully compared with two other important frequencies, i.e. 
the frequency of matter-radiation equality (be it v eq ) and the frequency of neutrino decoupling (which 
also coincides, in loose terms, with the Hubble radius at the onset of big bang nucleosynthesis). These 
two frequencies can then be written, respectively, as: 

„„ = ** = 1.281 x waggly h &™ r 1/2 Hz, (i.7) 

q 2vr V0.13267 V4.15x 10- 5 / v ; 

u hhn = 2.252 xloW-gg-V^f-^V ^ R ° , V /4 Hz^0.01nHz. (1.8) 

V10.75/ V MeVy V4.15x 10- 5 / V ; 



In Eqs. ( 11.71) and (11. 8p ^mo and Qro denote, respectively, the present critical fraction of matter and 
radiation with typical values drawn from the best fit to the WMAP 5-yr data alone and within the 
ACDM paradigm. In Eq. (11.81) g p denotes the effective number of relativistic degrees of freedom 
entering the total energy density of the plasma. While i/ eq is still close to the aHz, z/bbn is rather in 
the nHz range. 

The success of the CMB and BBN calculations implicitly demands that, after neutrino decoupling, 
the Universe was already dominated by radiation. If we assume that the radiation dominates right at 
the end of inflation, then the maximal frequency of the graviton spectrum can be computed and it is 
given by 

e y/ 4 / An \ 1/4 / ^Ro \ 1/4 



z/ max = 0.346 — ° ftU - GHz, (1.9) 

V0.01/ V2.41 x 10-V V4.15 x 10" 5 / ' V ; 

where An denotes the amplitude of the power spectrum of curvature perturbations evaluated at the 

pivot wavenumber k p . Between z/bbn and z/ max there are roughly 20 orders of magnitude in frequency. 

In the ACDM scenario the relic graviton spectrum has, in this range, always the same slope. 

1.2 CMB data and relic gravitons 

As already mentioned, CMB experiments are sensitive to long wavelength gravitons with typical fre- 
quencies of the order of v v ~ aHz (see also Eq. ( 11.11) ). The number of CMB parameters depends upon 
the specific model used to interpret (and fit) the data. The ACDM scenario probably contains the 
fewest number of parameters required to have a consistent fit of CMB data. 

The ACDM parameters can be inferred from various experiments and, among them, a central role 
is played by WMAP [51 El El El IS] (see also [ID| HIj H2j for first year data release and [El Ej for the 
third year data release) as well as other experiments (see, for instance, [15] in connection with the 5-yr 
WMAP data release). The TT, TE and, partially EE angular power spectrazl have been measured by 
the WMAP experiment. Other (i.e. non space-borne) experiments are now measuring polarization 
observables, in particular there are the 3-yr Dasi release [IE], the CAPMAP experiment [17], the recent 
QUAD data [19], [20], as well as various other experiments at different stages of development. The TT, 
TE and EE power spectra are customarily analyzed in the light of the minimal ACDM scenario but 
also other models are possible and they include, for instance, the addition of spatial curvature (i.e. 
the open-ACDM), more general parametrizations for the equation of state of dark-energy and so on 
and so forth. 



The combined analysis of the CMB data, of the large-scale structure data [221 E3] and of the 
supernova data [2U [25] can lead to quantitative upper limits on the possible contribution of the 
tensor modes to the initial conditions of the CMB temperature and polarization anisotropies. These 
upper limits can be phrased in terms of rx, i.e. the ratio between the power spectrum of tensor 



4 Following the custom the TT correlations will simply denote the angular power spectra of the temperature autocor- 
relations. The TE and the EE power spectra denote, respectively, the cross power spectrum between temperature and 
polarization and the polarization autocorrelations. 



fluctuations and the power spectrum of the scalar fluctuations evaluated at the pivot wavenumber 
k p = 0.002 Mpc™ 1 . In the minimal paradigm (i.e. the ACDM scenario) the tensor are not included in 
the fit. 

If the inflationary phase is driven by a single scalar degree of freedom and if the radiation dominance 
kicks in almost suddenly after inflation, n not only determines the tensor amplitude but also, thanks 
to the algebra obeyed by the slow-roll parameters, the slope of the tensor power spectrum, customarily 
denoted by n^. To lowest order in the slow- roll expansion, therefore, the tensor spectral index is slightly 
red and it is related to r^ (and to the slow-roll parameter) as n^ ~ — r^/8 ~ — 2e, where e = —H/H 2 
measures the rate of decrease of the Hubble parameter during the inflationary epoch Q Within the 
established set of conventions the scalar spectral index n s is given by n s = (1 — 6e + 2rj) and it depends 
not only upon e but also upon the second slow-roll parameter rj = M P V W /V (where V is the inflaton 
potential, V tW denotes the second derivative of the potential with respect to the inflaton field and 
¥p = l/V^G)- 

Depending upon the specific data used in the analysis, the upper limits on r T as well as the 
determination of the other cosmological parameters may change slightly. In Tab. 11.21 the upper 



Data 


r T 


n s 


^A 


^MO 


fc eq MpC 


WMAP5 alone 
WMAP5 + Acbar 
WMAP5± LSS + SN 
WMAP5± other CMB data 


<0.43 
<0.40 
<0.20 
<0.36 


0.986 ±0.22 

n qoco.oig 
u.yoj_ 0020 

0.968 ±0.015 

0.979 ±0.020 


77n+ - 033 
u -' 'U-0.032 

0.767 ±0.032 
0.725 ±0.015 
0.775 ±0.032 


230 - 032 
0.233 ±0.032 
0.275 ±0.015 
0.225 ± 0.032 


0.00936 
0.00944 
0.00999 
0.00922 



Table 1: The values of r^ are reported as they have been estimated in the absence of any running of 
the (scalar) spectral index. 



limits on rx are illustrated as they are determined from the combination of different data sets. For 
illustration the determined values of the scalar spectral index (i.e. n s ), of the dark energy and dark 
matter fractions (i.e., respectively, fi\ and ^mo), and of the typical wavenumber of equality k eq are 
also reported in the remaining columns. While different analyses can be performed, it is clear, by 
looking at Tab. I1.2l that the typical upper bounds on r T range between, say, 0.2 and 0.4. Slightly more 
stringent limits can also be obtained by adding supplementary assumptions. Within a conservative 
perspective, the tensor power spectra are, at least, ten times smaller than the power spectra of curvature 
perturbations. In the near future the Planck explorer satellite [26] might be able to set more direct 
limits on rx by measuring (hopefully) the BB angular power spectra |j. The E-mode power spectra and 
the B-mode power spectra arise as two orthogonal combinations of the Stokes parameters which are 
frame-dependent (i.e. Q and U). While the adiabatic mode leads naturally to the E-mode polarization, 
the only way of obtaining the B-mode (in the standard ACDM paradigm) is through the contribution 



5 The overdot will denote throughout the paper a derivation with respect to the cosmic time coordinate t while the 
prime will denote a derivation with respect to the conformal time coordinate r. 

forthcoming projects like Clover [27], Brain [25], Quiet [2j|] and Spider [3U] have polarization as specific target. 



of the tensor modes. Consequently, a detection o the BB angular power spectra would be equivalent, 
in the ACDM framework, to a first determination of r^. Having reviewed the essentials of CMB 
data and their connection with the relic graviton spectra, we can now move to higher frequencies and 
describe the status of the other devices which could shed light on the relic gravitons, i.e. the wide-band 
interferometers. 



1.3 Wide-band interferometers 

The wide-band interferometers operate in a window ranging from few Hz up to 10 kHz. The available 
interferometers are Ligo pQ, Virgo [2], Tama [3] and Geo jl]. The sensitivity of a given pair of wide- 
band detectors to a stochastic background of relic gravitons depends upon the relative orientation of 
the instruments. The wideness of the band (important for the correlation among different instruments) 
is not as large as 10 kHz but typically narrower and, in an optimistic perspective, it could range up to 
100 Hz. The putative frequency of wide-band detectors will therefore be indicated as z/lv, i-e. in loose 
terms, the Ligo/Virgo frequency. There are daring projects of wide-band detectors in space like the 
Lisa [31], the BBO [32J and the Decigo [33] projects. The common denominator of these three projects 
is that they are all space-borne missions and that they are all sensitive to frequencies smaller than 
the mHz. While wide-band interferometers are now operating and might even reach their advanced 
sensitivities along the incoming decade, the achievable sensitivities of space-borne interferometers are 
still on the edge of the achievable technologies. Since z/bbn < ^lv < ^max the wide-band interferometers 
are an ideal instrument to investigate the relic graviton spectrum in the unknown territory where there 
are neither direct nor indirect tests on the thermal history of the plasma. The problem is that, as it 
will be carefully shown, the spectral energy density of the relic gravitons produced within the ACDM 
model is quite minute and it is undetectable by interferometers even in their advanced version where 
the sensitivity is expected to improve by 5 or even 6 orders of magnitude in comparison with the present 
performances [3U [351 ES] (see also [37] and [38]). This impasse, as previously stressed, stems from 
the assumption that, right after inflation, the radiation-dominated evolution kicks in almost suddenly. 
At the moment, there are no evidences neither in favor of such a statement nor against it. The main 
theme of the present investigation will be to reverse this problem. It will be argued that wide-band 
detectors, in their advanced version, will be certainly able to test definite deviations from a simplistic 
thermal history of the plasma, i.e. the one stipulating that, after inflation, the radiation was suddenly 
dominating the evolution. 

1.4 Layout of the investigation 

In the present investigation the spectral energy density of the relic gravitons will be calculated first at 
small frequencies (compatible with the CMB observations) and then at higher frequencies (compatible 
with the operational window of wide-band interferometers). The latter calculation will be performed 
both n the case of the ACDM paradigm but also in those extensions which may lead to a large spectral 
energy density at the scale of the wide-band interferometers without violating any of the bounds 



stemming from CMB observations. 

The spectral energy density of the relic gravitons will be introduced in section [2j Gravity is 
inherently a non-Abelian gauge theory there are potential ambiguities in defining univocally an energy- 
momentum (pseudo)-tensor for the relic gravitons: it will be shown that different choices of the energy- 
momentum pseudo-tensor lead to the same spectral energy density of the relic gravitons. The punch 
line of section [2] will be that the spectral energy density can be more accurately performed with 
numerical methods rather than resorting, as often done, to semi-analytical estimates which amount to 
estimate first the power spectrum and then the spectral energy density. 

In section |3l using the numerical techniques introduced in section [2J the spectral energy density of 
the relic gravitons will be computed in the case of the standard ACDM paradigm which leads to nearly 
scale-invariant spectra. The signal arising in the context of the ACDM paradigm will be confronted 
with the sensitivity of wide-band interferometers. 

After showing the compatibility of the new methods with the results of the nearly scale-invariant 
spectra, possible scaling violations in the spectral energy density will be discussed in section HI While 
in the ACDM paradigm the spectral energy density of the relic gravitons is nearly scale invariant 
it is plausible to construct a class of models where the spectral energy density is fully compatible 
with the CMB and with the large-scale data at low frequencies while it is potentially detectable by 
wide-band interferometers. When we say that it is plausible this simply means that it is not forbidden 
by any of the current observational data. The proposed extensions of the ACDM paradigm have 
also a physical interpretation (see section H|) since they naturally arise when the thermal history of 
the Universe deviates, for sufficiently early times, from the usual assumptions of the ACDM scenario 
which stipulates that, right after inflation, the Universe suddenly becomes dominated by radiation. 

The minimal realization of the ideas pursued in section H] is scrutinized in section [5] and it is 
dubbed TACDM scenario (for tensor-ACDM). The TACDM paradigm consists of two supplementary 
parameters, i.e., in broad terms, a new pivotal frequency and a new spectral slope. The new frequency 
marks the onset of the high-frequency branch of the spectral energy density of the relic gravitons. In 
section [5] the TACDM paradigm is compatible with the current bounds stemming not only from CMB 
and large-scale structure. It will also be required that the big-bang nucleosynthesis (BBN) as well as 
pulsar timing constraints are satisfied and this will allow to spell out quantitatively the restrictions on 
the two supplementary parameters characterizing the TACDM scenario. 



2 Basic Equations 

The basic technical tools required to pursue the present analysis will be hereby summarized. The first 
part of the present section (i.e. subsection 12.11) contains an introduction to the basic terminology and 
a swift derivation of the main equations. Subsections 12.21 and 12.31 contain the details of our numerical 
approach whose results will also be illustrated in various physically relevant examples (see subsection 
12.41) and compared to the corresponding semi-analytical estimates (see subsection 12.51) . Finally, the 



exponential damping of the relic graviton spectrum will be numerically discussed in subsection 12.61 
As explained in the general layout of the investigation, all the considerations of the present section 
are bound to the ACDM paradigm so that the typical values of the cosmological parameters may be 
usefully drawn from Tab. 11.21 

2.1 Generalities 

In the ACDM paradigm the geometry is conformally flat (see Eq. (11.21) ) and the corresponding tensor 
fluctuations are defined as 

Si 1] 9ij = -a\r)h lv 8?g« = ^, 5^ g* = -^, (2.1) 

where h\ = 8^ = 0. The second order action obeyed by hy can be written as 

1 
-p 

where £p is defined as Q 



S GW = ^f I d 4 x^g^d,h tJ d u hv, (2.2) 



t = 8ttG = i = |1. (2.3) 

P M P Ml K > 

Equation ( 12.21) is effectively equivalent to the sum of the actions of two (scalar) degrees of freedom 
minimally coupled to the background geometry. To derive Eq. (12. 2p the Einstein-Hilbert action must 
be perturbed to second order in the amplitude hij, i.e. 

6 ^ S = -^h [d 4 x[S? ) V^9R+V^5i 2) R + 5[ 1) ^g5i 1) R]. (2.4) 

lo7rG J 

To evaluate Eq. (12.41) in explicit terms it is necessary to compute the Ricci tensors both to first and 
second order, i.e. 

^ ] R i0 = ^[tit j + 2Hti ij -V% j ] + (H' + 2H 2 )h ij , (2.5) 



f( 2 ) E>__ — _UI h^' _ _h..hij' _L _MfV72j 

4 
1 



S[ 2) R 00 = -h! l3 W' - -h tJ h^' + -tiWhij, (2.6) 



4 v 2 J 2 

f(2)D.. _ l hki\ 



^t R%j — tJ 11 [dkdehij — dkdjhn — d k dih 3 



I T-/ 

-dj[h ke (d e h ik - d k h u - dih H )} - —i 



ft i/i/ l-,£u' Uk' "ul Uk' 'u Uk 1 J,/ uk 1 u 

—tijtiti + —hihfj - -hj n ik - —hj n ik - -fy fyy - — ^ /i fcj - 

i[9^i + ^M - ^ fc ][9//iJ + d 3 h1 - 9 fc M- (2-7) 



7 Some authors include a V87T is the definition of the reduced Planck mass (what we call Mp = £p ). This convention 
is per se harmless, however, it may be confusing in practice. In the present script the conventions expressed by Eq. 
will always be carefully followed. 



The Ricci scalar is zero to first order in the tensor fluctuations, i.e. <5 t R = 0. This is due to the 

(2) 

traceless nature of these fluctuations. To second-order, however, 61 R ^ and its form is: 



W R = i^l^^' + ^^^ + I^v 2 ^-^^^^ 



^[dihidth* - d t hid k h u + d k h a d e h k - ^hikd'hl + d% k d k h u \y (2.8) 

Using the results of Eqs. (I2.5I) - (I2.8I) into Eq. ( 12.41) the second-order action for the tensor modes of 
Eq. (I2.2p can be obtained by getting rid of a number of total derivatives. 

According to Eq. (12.11) . hij carries two degrees of freedom associated with the two polarizations of 
the graviton in a Friedmann- Robert son- Walker (FRW) space-time. Defining as k = k/\k\ the direction 
along which a given tensor mode propagates, the two polarizations can be defined as 

e®-(fc) = (rhirhj - h^), e|(fc) = (rhifij + rhjhi), (2.9) 

where rh and n are two mutually orthogonal unit vectors which are also orthogonal to k (i.e. m ■ n = 
n ■ k — m • k — 0). During the early stages of the ACDM model (i.e. during the inflationary phase) 
hij(x,r) can be expanded in terms of the appropriate creation and annihilation operators as: 

M*, r) = ^j- 2 E / d 3 k ef{k) [F M (T)o SA e-** + Fl x {r)a\/ k % (2.10) 

where the index A counts the two polarizations, i.e. A = ®, ©; k denotes the wavenumber and F k \(r) 
is the (complex) mode function obeying 

Kx = G ^ ( 2 - n ) 

G' KX = -2HG kA -k 2 F k>x , (2.12) 

In Eq. ( 12.101) [a kX ,al x ,] = 6^ 3 \k — p)6 XX i. The initial state |0) (annihilated by a^ A ) minimizes the 
tensor Hamiltonian when all the wavelengths of the field are shorter than the event horizon at the onset 
of the inflationary evolution. The main observables which are used to characterize the relic graviton 
background are the two-point function evaluated at equal times and the spectral energy density in 
critical units. The two-point function is defined as 

-s -s /*oo sin k'T' 

(0\h ij (x,T)h ij (it,'rM = J dlnkV T {k,r)-^—, r = \x - y\, (2.13) 

V T (k,r) = ^#|F fc (r)| 2 . (2.14) 

The quantity V^{k, r) is, by definition, the tensor power spectrum. Equations (I2.13l) - fl2.14l) can be 
derived by recalling the following pair of relations: 

ejfejf = 26 xy , F k>(B (r) = F k>9 (r) = F k (r). (2.15) 

8 



Out of Eqs. (I2.13p - (l2.14p it is sometimes practical to introduce the spectral amplitude Sh{y, t)j namely, 

(0\h ij (x,T)h ij (x,r)\0)= V T {k,r)dlnk = 4 / vS h (u,r)d\nu, (2.16) 

where fe = 2iru. By definition, pq^(x,t) = (0|Tq (x, r)|0) where |0) is, again, the state annihilated by 
afcA (see also Eq. (12.101) ) and where T^ is the energy-momentum (pseudo)-tensor of the relic gravitons. 
The spectral energy density of the relic gravitons in critical units can then be computed from the 
expectation value of Tq lj 

n G w(k,T) = —^-, p GW = (0|T |0>. (2.17) 

p crit dink 

where p cr i t = 3H 2 /£p is the critical energy density. In FRW space-times the energy-momentum pseudo- 
tensor of the relic gravitons can be assigned in manners which are conceptually different but physically 
complementary. This will be one of the topics discussed in subsections 12.21 and 12.31 

2.2 Transfer function for the amplitude 

To connect the early moment of the normalization of the relic gravitons to the moment where the 
tensor modes of the geometry reenter the Hubble radius and affect, in principle, terrestrial detectors 
the customary approach is to solve for the tensor mode function and to define the so-called amplitude 
transfer function. From the amplitude transfer function the spectral energy density (see e.g. Eq. 
(I2.17P ) can be computed. What we propose here is to do the opposite, i.e. to compute, numerically 
and in one shot, the transfer function for the spectral energy density. To show the equivalence (but 
also the inherent differences) of the two approaches, the present subsection will be concerned with 
the transfer function for the amplitude. In the following subsection (i.e. subsection I2.3P the transfer 
function of the spectral energy density will be more specifically discussed. 

During the inflationary phase, the tensor power spectrum can be easily computed by solving Eqs. 
(12. lip and (I2.12p in the slow-roll approximation 

F*(t) = -"v=ktHW(-kr), M = JZ e *MW, u = -^-. (2.18) 

a\T)\J2k V z Z{L — e) 

where HJP(z) = J v {z) +iY v (z) is the Hankel function of first kind [391 HO] and where e = —H/H 2 . To 
obtain the result of Eq. (12.181) from Eqs. (12.111) and (I2.12p it is useful to bear in mind the following 
pair of identities 

H 2 + H 2 = a 2 H 2 (2-e), aH = -^ — -. (2.19) 

r(l -e) 

The second equality in Eq. (12.191) can be simply deduced (after integration by parts) from the relation 
between cosmic and conformal times, i.e. a{r)dr = dt. Physically Eqs. (12.181) and (12.191) hold under 
the approximation that, at early times, the background geometry is of quasi-de Sitter type. 



3 Thc natural logarithms will be denoted by In while the common logarithms will be denoted by log. 



By substituting Eq. (12.181) into Eq. (12.141) the standard expression of the tensor power spectrum 
can be obtained. When the relevant modes exited the Hubble radius during inflation: 

V T (k, r) = ilH^(u)(l - C) 2 "- 1 (^)^, v = \ + e + 0(e 2 ), (2.20) 

where the small argument limit of the Hankel functions has been taken and where the slow-roll ap- 
proximation has been enforced, i.e. in formulas: 

k k H , , N 

x - kT ^H = ^H <<1 ' e = -JP <L (2 ' 21) 

The two approximations introduced in Eq. (12.211) will be often employed and, therefore, it is appro- 
priate to spell out clearly their physical meaning. The first relation of Eq. (I2.2ip implies that kr < 1 
this means that the wave-numbers are, in practice, smaller than the Hubble rate. Conversely the 
corresponding wavelengths will be larger than the Hubble radius. The chain of equalities appearing 
in Eq. (12.211) can be easily understood since, by definition, 7i = aH and H is the Hubble rate. In 
this long wavelength limit, as we shall see, the evolution of the tensor modes can be derived in semi- 
analytical terms and it corresponds to a tensor mode function F^{r) which is approximately constant. 
The latter statement holds if the geometry is of quasi-de Sitter type. The latter condition is verified 
if the second relation of Eq. (I2.2ip . stipulating e < 1, holds. The latter conditions is also dubbed 
slow-roll approximation and it allows to simplify the tensor power spectrum even further: 

— 2 / V \ 128/ V \ 

^-wK^f) -T-(Xfi) • (2 ' 22) 

OK \MpSk~aH 3 \M^/k~aH 

The spectral index defined from Eq. d 2 . 2 2 [) is nothing but 

" T = l^ = -T^ = - 2£ + ^- (2 ' 23 > 

where the second equality can be derived with the standard rules of the slow-roll algebra. The spectral 
amplitude and slope are then parametrized, for practical purposes, as 

p T (fc)=^ T f_J , k p = 0.002 Mpc- 1 , (2.24) 

where, by definition, At is the amplitude of the tensor power spectrum evaluated at the pivot scale 
k p . The pivot wavenumber of Eq. (I2.24p is simply related to the pivot frequency defined in Eq. (11.11) 
as Up = fc p /(27r). Bearing in mind that the power spectrum of curvature perturbations is given, in 
single field inflationary models, as 

p « (fc >=3taj) M =Ms-) • (2 - 25) 

the ratio between the tensor and the scalar power spectra is simply given by 

7*r(fc) A T 

rT = rrii = ~r~ = 16e ' 2 - 26 
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Equation (I2.26J) implies, recalling Eq. fj2.23|) . that r^ = — 8tit- In Tab. 11.21 the values of rx have been 
reported as they can be estimated in few different analyses of the cosmological data sets. 

Equation (12.241) correctly parametrizes the spectrum only when the relevant wavelengths are larger 
than the Hubble radius before matter-radiation equality. To transfer the spectrum inside the Hubble 
radius the procedure is to integrate numerically Eqs. f l2.lip - fl2.12l) (as well as Eqs. fll.4j> — (11.6D ) across 
the relevant transitions of the background geometry. While the geometry passes from inflation to 
radiation Eq. (12.241) implies that the tensor mode function is constant while the relevant wavelengths 
are larger than the Hubble radius: 



F t{ r) = A t + B t j^-r ±«1, |^ = 3 |1jP tW . (2.27) 



The term proportional to B k in Eq. (12.271) leads to a decaying mode and Fk(r) is therefore determined, 
for \kr\ <C 1, by the first term whose squared modulus coincides with the spectrum computed in Eq. 
(12.22p and parametrized as in Eq. (12.241) . The evolution of the background (i.e. Eqs. (jl.4]l — (J1.6D ) and 
of the tensor mode functions (i.e. Eqs. fl2.lip - fl2.12p ) should therefore be solved across the radiation 
matter transition and the usual approach is to compute the transfer function for the amplitude 
i.e. 



T h (k) 



\F k (r W ) ( ^ 2g) 



\ (ln(r)| 2 ) 



In Eq. ( 12.28ft . Fk{T~) denotes the approximate form of the mode function (holding during the matter- 
dominated phase); F k {j) denotes, instead, the solution obtained by fully numerical methods. As the 
wavelengths become shorter than the Hubble radius, F k {r) oscillates. Consequently, To get T h (k) the 
oscillations must be carefully averaged and this is the meaning of the averages appearing in Eq. f)2.28f) . 
Hence, the calculation of Th(k) requires a careful matching over the phases between the numerical and 
the approximate (but analytical) solution. Consider, indeed, one of the most important applications 
of the previous results, i.e. the radiation-matter transition. After matter-radiation equality, the scale 
factor is going, approximately, as a{r) ~ r 2 and, therefore, the (approximate) solution of Eqs. (I2.1ip - 

(I2.12p is given by 

— . . 3?i(Arr) . ., . sinfcr cosfcr 

*■»« = ^r 1 ^. #) = j^ - w (2.») 

which is constant for kr < 1. In Fig. [I] the result of the numerical integration is reported in terms of 
T%(k/k cq ). In Fig. [I] the fit to the numerical points is also reported and it can be parametrized as: 



T h (k/k 



cq/ 



l + Cl (A)+ 6l (M 2 . (2 .3 0) 



By applying the standard tools of the regression analysis c\ and b\ can be determined as c\ = 1.260 
and b\ = 2.683. The latter result agrees with the findings of jH] who obtain c\ = 1.34 and b\ = 2.50. 
The value of k eq can be obtained directly from the experimental data (see, for instance, last column 
of Tab. 11.21 implying k eq ~ (9(0.009) Mpc~ ). For instance, the WMAP 5-yr data combined with the 
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5J£ Numerical integration 
^— ^— Fit to the numerical data 



C\J JZ 



CO 

o 



» (: | c) | c^)^ | ti | ;i | ;- | w | w | :: | i: | H* 




iog(k/k eq ) 



Figure 1: The starred points represent the numerical values of the amplitude transfer function of the 
amplitude across the matter-radiation transition. The logarithm (to base 10) is reported on both axes. 
The full line represents the numerical fit. 



supernova data and with the large-scale structure data would give k eq = 
turns out that a rather good analytical estimate of k eq can be presented as 



n nnQQQ+ 000028 
u.uuyyy_ 00027 



Mpc" 



It 



h 



eq 



0.0082879 



^0^M0 

0.1326 



/U.15x 10 



-1/2 



Mpc" 



(2.31) 



where the typical value selected for h^fl^Q is given by the sum of the photon component (i.e. /ig^o = 
2.47 x 10~ 5 ) and of the neutrino component (i.e. ^0^70 — 1-68 x 10~ 5 ): the neutrinos, consistently 
with the ACDM paradigm, are taken to be massless and their (present) kinetic temperature is just a 
factor (4/11) 1 / 3 smaller than the (present) photon temperature. From Eq. (I2.3ip it is straightforward 
to estimate the equality frequency of Eq. (jl.7p . 

The analytical estimate stems from the observation that the exact solution of Eqs. (jl.4p - (jl.6p for 
the matter-radiation transition can be given as a(r) = a eq [y 2 + 2y] where y = t/t\. The time-scale 
T\ = r eq (\/2 + 1) is related to the equality time r eq which can be estimated as 



2(^-1)^ 



R0 



' eq 



Hn 



n 



120.658 



M0 



0.1326/ 



hltt 



O^RO 



4.15 x 10 



-5 



1/2 



Mpc. 



;2.32) 



In the case of the WMAP 5-yr data combined with the supernova and large-scale structure data 
h 2 Q M0 = 0.1368™° 3 o 8 37 _ Consequently, Eqs. (J22HD, jOg) and fl23UD imply that the spectrum of the 
tensor modes is given, at the present time, as 



^x(A:,r ) 



(kr ) 



T 2 h (k/k cq )V T (k). 



(2.33) 



Within the standard approach, Eq. (I2.33P is customarily connected to the spectral energy density of 
the relic gravitons. It will now be shown that the spectral energy density of the relic gravitons can be 
directly assessed, by numerical means, without resorting to Eq. fj2.33|) . 
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2.3 Spectral energy density 

Having presented the standard derivation of the amplitude transfer function we will now focus the 
attention on the transfer function of the spectral energy density. The latter approach leads more 
directly to the estimate of the present value of the spectral energy density. Before discussing in some 
detail the numerics it is appropriate to recall the construction of an energy-momentum pseudo-tensor 
by following the same approach which has been proven to be successful in flat space-time [32] ■ in a 
conformally flat geometry of the type introduced in Eq. (jl.2p . the energy momentum pseudo-tensor can 
be derived by following two complementary strategies. The first one is to take the energy-momentum 
tensor associated with the action of Eq. ( 12. 2D . Since each polarization of the graviton in a FRW space- 
time obeys the evolution equation of a minimally coupled scalar field, it is legitimate to establish that 
the energy- momentum pseudo-tensor is just given by the energy- momentum tensor of a pair of scalar 
degrees of freedom minimally coupled to the geometry. By formally taking the functional derivative 
of Eq. (12.21) with respect to g^ u , T^ becomes 



rpv 

» ~ Ml 



1 



d^Ptf - -Stftdahi&hV = -s E 



2 



1 



U 



d»h ix) d»h^ - \r p d a h {x) d p h w 5; 



(2.34) 



p 



where the second equality follows from the first by using that h^ = Y,x ^{\) e ij an d that e\j e\j = 25 W>. 
This perspective was adopted and developed, for the first time, in [4"3"t 133] by Ford and Parker. A 
complementary approach is to use the energy-momentum pseudo-tensor defined from the second-order 
fluctuations of the Einstein tensor: 

V = ~^?% Q v , = Rl - \^R- (2.35) 

(2) 

where <5 t denotes the second-order tensor fluctuation of the corresponding quantity. The latter ap- 
proach is more directly related to the well known flat space-time procedure [32] • The approach ex- 
pressed by Eq. ( 12.351) has been described in [33, I3S] and has been reprised, in a related context, by 
the authors of Refs. [33 I3H] mainly in connection with conventional inflationary models where the 
Universe is always expanding. 

According to Eq. (12.341) the energy density is given by p G ^ = (0|Tq |0) where |0) is the state 
annihilated by the creation and destruction operators introduced in Eq. (I2.10p : 

&(r) = ^ J dlnk^- 2 {\g k (r)\ 2 + (k 2 + H 2 )\f k (r)\ 2 - H[r k (r)g k (r) + f k (r)g* k (r)]], (2.36) 

where fkij) = F k {T)a{r) and guij) = f'k( T ) have been introduced. The superscript appearing in Eqs. 
(I2.36P reminds that the energy density refers to the first choice of the energy-momentum tensor given 
in Eq. ( 12. 34ft . According to Eqs. ( 12.111) and ( 12.121) the tensor mode functions f k and g k obey 

f' k = 9k, g' k = -[k 2 - (H 2 + H')]f k . (2.37) 

By adopting the approach expressed by Eq. (j2.35|) . the energy density of the relic gravitons p G ^ = 
(0|T °|0) become: 

&(r) = J dlnk ^l{\9^)\ 2 + i k " ~ 7n 2 )\f k (r)\ 2 + 3H[r k (r)g k (r) + /*(r)^(r)]}. (2.38) 
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To pass from Eq. (I2.35J) to Eq. (J2.38P the simplest procedure is to obtain the second-order fluctuation 
of the Einstein tensor, i.e. 5^Q U . This calculation can be easily carried on by using the results of Eqs. 
(I2.5I) - (I2.8I) . Furthermore, it should be appreciated that the two energy- momentum pseudo-tensors 
lead also to different pressures and this observation has an impact on the back- reaction problems |49j . 
From Eqs. ( 12.361) — ( and (12.381) the corresponding critical fractions are: 



^Gw(^ T ) 



1 dp, 



(i) 

GW 



"Gw(^ r ) 



1 dp 



(2) 
GW 



p crit d In k ' ' p crit d In k 

If k/Ti > 1, then /^(t) will be, in the first approximation, plane waves and ^(r) 
two versions of Hgw(^ t ) will be given by: 



(2.39) 
=t^/fe(' r ) and the 



»(i) 



^Gw(^ r ) 



"Gw("')' r ) 



kH\ 
37T 2 a 2 H 2 

kH\ 
37v 2 a 2 H 2 



1 + 



2k 2 

m 2 

2k 2 



\f k (r)\ 2 - 
\h(r)\ 2 



k 2 



12H 2 
k 2 



12H 



Vr(k,r) 

jV T (k,T 



2f_- 

2k 2 : 

7H 2 

' 2k 2 



(2.40) 
(2.41) 



where, the second equality follows from the first by recalling that \fk{ T )\ 2 — 7T 2 a 2 VT(k,r)/(A£pk 3 ) 
Equations (12.401) and (12.411) coincide (up to corrections OiTL 2 /k 2 )). This means, physically, that the 
energy density of the relic gravitons is effectively the same no matter what choice of the energy- 
momentum pseudo-tensor is adopted but provided the wavelengths of the gravitons are all inside (i.e. 
shorter than) the Hubble radius. When the given wavelengths are larger than the Hubble radius (i.e. 
fer < 1), jit = Hfk and the corresponding expressions for the spectral energy densities can be easily 
obtained. In summary, for modes which are inside the Hubble radius the energy density of the relic 



h =0.719, h„n,.n =0.1326, f!,„ + Q. =1 

MO MO A 



CO 

o 
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logx 
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Figure 2: The functions given in Eqs. ( 12.451) and ( 12.461) are numerically computed (plot at the left) 
for different values of k and in the case of the radiation-matter transition. In the plot at the right the 
transfer function for the energy density is illustrated. 



gravitons can be expressed in terms of the power spectrum as 

pow(T ) =^/«*ln*^|A|» = ^/dh**»p r (*,r) 
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(2.42) 



and the critical fraction of relic gravitons at a given time as: 

^gw(^) = — p% = -*V?(k,r). (2.43) 

p crit a ink 12H z a^ 

Specific examples of the numerical calculation of the spectral energy density will be given in the 
following subsection (i.e. subsection 12.41) . The first example will be the one of the radiation-matter 
transition the second example will be the one of the transition between a stiff phase and the radiation- 
dominated phase. 



2.4 Transfer function for the spectral energy density: examples 

The idea pursued in the present subsection is to use, as pivot quantity for the numerical integration, 
not the power spectrum V(k, r) but rather the energy density itself. The evolution equations of the 
background geometry (i.e. Eqs. fll.4l) - (ll.6p ) and of the tensor mode functions (i.e. Eq. (12.37P ) will be 
solved simultaneously and the energy density computed in one shot. This program will be illustrated 
in two simple examples, i.e. the radiation-matter transition and the case of a stiff background. It is 
useful to point out that Eq. (12.211) suggests that an appropriate variable for the numerical calculation 
is exactly x = kr whose definition we now repeat: 

x = kr = k( — J, k = — . (2.44) 

\ 'eq ' ""eq 

It is both practical and physically sound to adopt x and k as pivotal variables for the numerical 
integration around the radiation-matter transition. Indeed x is a smooth variable which interpolates 
between the sub-Hubble regime (where Eqs. (I2.40I) - (I2.41I) are valid and the super-Hubble regime 
where x > 1. The result of the numerical calculation are reported in Fig. [2] in terms of Ay(K,x) and 
in terms of the transfer function of the energy density (denoted by T p (k)). The quantities AW(k, x) 
(and, analogously, Ay(K, x)) are nothing but 

A?> (*,T) = {k(r)| 2 + (fc 2 + H 2 )|/ fc (r)| 2 -W[/;(r) ft (r) + /,(rM(r)]}, (2.45) 

A?>(*,t) = {k(r)| 2 + (fc 2 -7H 2 )|/ fe (r)| 2 + 3H[/*(r)^(r) + / fc (rK(r)]}, (2.46) 

Equations (I2.45J) and (12.461) are simply related to the spectral energy densities in critical units, i.e. 

«SU*,t) = - f, AW(fc,r), n^(*,r) = - f, Af\k,r). (2.47) 

27r J a 4 p crit F 27r i a 4 p crit p 

As a function of x and k A^' 2 \k, x) reaches a constant value when the relevant modes are evaluated 
deep inside the Hubble radius. The energy transfer function which is then defined as: 

lim A£' 2 \k,x) = T*(k)A£' 2 \k,x } ), x, < 1. (2.48) 

The specific form of the energy- momentum tensor is immaterial for the determination of T 2 (k): 

15 
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Figure 3: The different definitions of energy-momentum pseudo-tensor (i.e. Eqs. (I2.45I) and (12.460 ) 
are compared in the determination of the asymptotic value of the energy transfer function. 



different forms of the energy-momentum tensor of the relic gravitons will lead to the same result. This 
occurrence can be appreciated from Fig. [3] where A^' 2 \k, r) has been reported for k = 10~ 2 (plot at 
the left) and for k = 10~ 4 (plot at the right). The dashed and the dot-dashed curves (in both plots) 
correspond, respectively, to A^'(k,x) and to A^(k, x). The full line, in both plots, corresponds to 
the combination 

k 2 \fk(r)\ 2 + \9k(r)\ 2 = H\c + (k)\ 2 + |c_(fc)| 2 ), (2.49) 

where c±(k) are the so-called mixing coefficients which parametrize, at a given time, the solution for 
the tensor mode functions when the relevant wavelengths are all inside the Hubble radius, i.e. 



hi?) 



c + (k)e- tkT + c_(fc)e' 



Ikr 



9k{T) 



Ik 

M/ 2 



c + {k)e 



-ikr 



c_(k)e 



ikr 



(2.50) 



where f k (r) and g k (r) are the solutions to leading order in the limit kr ^> 1. From Eq. (j2.50p . c±(k) 
are given by 



c + {k) 



/-*%r\>t /-* XrZi 



[kf k (r)+ig k (r)}, c^k) 



[kfk(T)-ig k (r)], 



(2-51) 



'2k L v ' V2k 

Using Eqs. (12.50H — ( l2.51[) . Eqs. ( I2.45l) - (12.46l) can be directly assessed in the limit x = kr 3> 1 with the 
result that 

A^(K,x i ) = Af\K,x l )=K(\c + (K)\ 2 + \c4K)\ 2 )+0^, (2.52) 

which proofs that the oscillating contributions are suppressed as a;^ 1 for Xf ~^> 1. 

To get to the results illustrated in Figs. [2] and [3] the evolution equations of the mode functions 
have been integrated by setting initial conditions deep outside the Hubble radius (i.e. x = kr <C 1), 
by following the corresponding quantities through the Hubble crossing (i.e. x ~ 1) and then, finally, 
deep inside the Hubble radius (i.e. i>1). The initial value of the integration variable x has been 
chosen to be x\ = 10~ 5 . The integration of the mode functions is most easily performed in terms of 
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appropriately rescaled variables and since these rescalings are rather obvious, the relevant details will 
be omitted. 

In the plot at the right of Fig. [2j the fit to the energy transfer function is reported with the full 
(thin) line on top of the diamonds defining the numerical points. The analytical form of the fit can 
then be written as: 



2 



T p (k/k eq ) = Jl + c 2 (-jl\ + b 2 (-jA , c 2 = 0.5238, b 2 = 0.3537. (2.53) 

Equation (12.531) permits the accurate evaluation of the spectral energy density of relic gravitons, for 
instance, in the minimal version of the ACDM paradigm. 

Yet another relevant physical situation for the present considerations is the one where the back- 
ground geometry, after inflation, transits from a stiff epoch to the ordinary radiation-dominated epoch. 
In the primeval plasma, stiff phases can arise: this idea goes back to the pioneering suggestions of 
Zeldovich [50] in connection with the entropy problem. The approach of Zeldovich was revisited in 
[5Tj |52| |53| [53] by supposing that the stiff phase would take place after the inflationary phase with 
the main purpose of identifying a potential source of high-frequency gravitons which could even be 
interesting for the LIGO/VIRGO detectors in one of their advanced versions. 

At the end of inflation, in a model-independent approach, it is plausible to think that the onset of 
the radiation-dominance could be delayed. This may happen, in concrete models, for various reasons. 
One possibility could be that the inflaton field does not decay but rather changes its dynamical nature 
by acting as quintessence field [59] (see also [60]). In this kind of situations the geometry passes from 
a stiff phase where 

Mr) = ~>i (2-54) 

Pt 3 

o p[ wi ldln(iu t + l) 1 , nrr , 

<%t(r) = -T-^t- ^77 \ u =w t -- f. ->-, 2.55 

p\ 3n(wt + l) 3 dma 3 

to a radiation-dominated phase where c st = l/\/3. According to Eqs. (12.541) and (12.551) . c^ t = w t iff 
the (total) barotropic index is constant in time. In the limiting case w% — 1 = c[? t and the speed of 
sound coincides with the speed of light. As argued in [57], barotropic indices w% > 1 would not be 
compatible with causality (see, however, [58]). As in the case of the matter-radiation transition the 
transfer function only depends upon k which is defined, this time, as k = k/k s , where k s = r" 1 and r s 
parametrizes the transition time. A simple analytical form of the transition regime is given by 



a(y) = a s Jy 2 + 2y, y = — , r s = —-J-, (2.56) 

v t s aiHi y p Ri 

where, by definition, p si = p s (ri) and pri = Pr(t ; ). Equation ( 12.561) is a solution of Eqs. fll.4p - fll.6D 
when the radiation is present together with a stiff component which has, in the case of Eq. ( I2.56P a 
sound speed which equals the speed of light. In the limit y — *■ the scale factor expands as a(y) = \f2K) 
while, in the opposite limit, a(y) ~ y. In Fig. [4] (plot at the left) A p (k,x) is illustrated for different 
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values of k. We shall not dwell here (again) about the possible different forms of the energy momentum 
pseudo-tensor. The bottom line will always be that, provided the energy density is evaluated deep 
inside the Hubble radius the different approaches to the energy density of the relic gravitons give the 
same result. From the numerical points reported in Fig. H] (plot at the right) the semi-analytical form 
of the transfer function becomes, this time, 

T 2 p (k/k s ) = 1.0 + 0.204 (j\ - 0.980 (£\ * + 3.389^) - 0.067 (y) In 2 (k/k B ), (2.57) 

where k s = r^ 1 . The value of k s can either be computed in an explicit modeo or it can be left as a free 
parameter. In section [3] both strategies will be explored by privileging, however, a model- independent 
approach. Taking into account that the energy density of the inflaton will be exactly p S i ~ H?M P , the 
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Figure 4: The transition between the stiff phase and the radiation phase is illustrated. The energy 
transfer function increases with the frequency while the opposite is true for the radiation-matter 
transition (see Fig. [2]). 

value of k s (as well as the duration of the stiff phase) will be determined, grossly speaking, by H[/Mp. 

2.5 Analytical estimates of the mixing coefficients 

To obtain a fit of the transfer function for the spectral energy density it is useful to be aware of the 
analytical results which should always be reproduced by the numerical analysis when k is sufficiently 
larger than 1. This is the purpose of the present subsection where it will be shown that the semi- 
analytical results are consistent with the numerical evaluations which are, however, intrinsically more 
accurate. 

Consider the transition from a generic accelerated phase to a decelerated stage of expansion. In 
this situation, by naming the transition point — r l5 the continuous and differentiable form of the scale 



In the context of quintessential inflation J5§] (see also 



J) m - Ht 



factors can be written as: 
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T < -Ti, 



a \n 



t > -Ti, 



(2.58) 
(2.59) 



where the scale factors are continuous and differentiable at the transition point which has been gener- 
ically indicated as T\. The pump fields of the tensor mode functions turn out to be: 



of _ f3({3 + 1) 
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iTl 



;2.60) 



The solution of Eq. (12.371) can then be written as: 
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/s(r) 



/2fc 

v 7 ^ 



^(- 



.r 



r < -Ti, 
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(2.61) 



where y = kr + feri ( 1 + f ) and where 
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;2.62) 



The continuity of the tensor mode functions at the transition point [i.e. fi(-Ti) = / s (— T\) and 
9\{~ T i) = 9s{— T i)} implies that the mixing coefficients are given by: 



c + {k) 



c-{k) 



in 



itt(u+X)/2 



tap 



{[(3{2X + 1) + a(2u + l)]Hi 1 \x 1 )H^(y 1 



r(l) 



2ax 1 [H^(y 1 )Hl 1 l 1 (x 1 ) + H^(x 1 )H^ 1 (y l )}} 



(i) 



in 



Jn(v-\)/2 



'af3 



{[(3{2\ + 1) + a (2u + l)]Hi l \x 1 )H^(y 1 



r(l) 



2ao; 1 [tf} 2 V)^ 1 + ) 1 (* 1 ) + ^(n)^^)]} 
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(2.63) 



where, according to the notations previously established, y\ = y{—T\) = {a/ (5)x\. The case a = (3 — 1 
corresponds to a transition from the inflationary phase to a radiation-dominated phase. In this case 
we do know which are the mixing coefficients. The previous expressions give us: 



c.(k) 



_2ixi 



2a;?' 



c+(fc)= ( i - 2 4 + ^" 



(2.64) 



which clearly agree with previous results [621 [63]. In the case of Eq. (12.641) |c + (fc)| 2 — |c_(/c)| 2 = 1 and 
k A \c-(k)\ 2 is exactly scale- invariant. Another interesting situation is the one of the transition from 
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inflation to stiff, i.e. /3 = I, a = 1/2, y\ = X\j2 which leads to a logarithmic enhancement at small 
wavenumbers [5TL 152] . In this situation the mixing coefficients can be written as: 

+ (i + x 1 )[ix 1 H[ 2 \x 1 /2)-3iH^\x 1 /2)]y (2.65) 

c + (k) = ^^^e^{x 1 4 1 \x 1 /2) + i(i + x 1 )Hi 1 \x 1 /2)y (2.66) 

The above result can be expanded in for X\ <C 1 and the result is: 



3/2 >, 



c+(k) = — —i - + ^xl [(0.131 + 0.338 i)- 0.149 (1 - i) lnxi] + C(xi /Z ), (2.67) 

X\* 



a_(*) 



(7.031 - 1.723*) - 16.68 (1 + i) \a.x x 



3/2 
X i 



+ V%I [(-0.621 + 0.265 i) + 0.282 (1 + i) In x x \ + O^' 2 ). (2.68 



The logarithms arising in Eqs. (12.671) and (I2.68P explain why, in Eq. (I2.57p . the transfer function of 
the spectral energy density contains logarithms. In spite of the fact that semi-analytical estimates 
can pin down the slope of the transfer functions in different intervals, they are insufficient for a 
faithful account of more realistic situations where the slow-roll corrections are relevant and when other 
dissipative effects (such as neutrino fee streaming) are taken into account. 

2.6 Exponential damping of the mixing coefficients 

In a model-independent perspective, it can be argued that the relic gravitons are also characterized by 
a maximal frequency which is related to the modes which are maximally amplified. Let us consider, 
for instance, the case of the ACDM paradigm where the inflationary phase is almost suddenly followed 
by the radiation-dominated phase. By denoting the transition time as r i5 it is plausible to think that 
all the modes of the field such that k > a^Hi ~ r^ 1 are exponentially suppressed [6U [65]. For the 
modes kv x > 1, the pumping action of the gravitational field is practically absent. The wavenumber 
km&x (which is related to the maximal frequency introduced in Eq. (11. 9p ) is the maximally amplified 
wavenumber which can be determined by requiring k ~ rf 1 : 



fc max = 3.5661 x 10 22 — ^^7 ° r ) Mpc-\ (2.69) 

V0.01/ V2.41xl0- 9 / V4.15xl0- 5 / P ' K J 

where the typical values of the slow-roll parameter have been derived by taking into account that, in 
the absence of running of the tensor spectral index, r^ = 16e; since, according to the WMAP 5-yr 
data alone, tt < 0.43, e < 0.01. Note that ^ max = 2nk maiX = 117.45 x (ireAn) 1 ^ GHz for the same 
typical values of the ACDM parameters. For phenomenological purposes it can be also interesting to 
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Figure 5: The time evolution of the mixing coefficients is reported at the left (on the horizontal axis 
the scale is linear). The exponential decay of the mixing coefficients is illustrated in the plot at the 
right. 

know what kind of exponential suppression we can expect. From the analysis of various transitions it 
emerges that the mixing coefficients for k > k max (or v > f max ) will satisfy 

-2/3 tt-^ 

P ^ Km a. 



|c + (fc)| 2 -|c_ 



M*)|' + |c_(*)|' 



+ 1. 



(2.70) 



From Eq. (I2.70p we can easily argue that, for k > k max , \c + (k)\ — > 1 and |c_(/e)| ~ 2 -1 / 2 exp [—fik/k max ]. 
The point is then to estimate the value of /3 which depends on the nature of the transition regime. 
Typically, however, (3 > 2 for sufficiently smooth transitions. To justify this statement it is interesting 
to consider the following toy model where the scale factor interpolates between a quasi-de Sitter phase 
and a radiation-dominated phase: 

a(r) = a i [r + v /r 2 + 7i 2 ]. (2.71) 

For r — > — oo (i.e. r <^ — r ; ) , a{r) ~ —cl-Jt and the quasi de-Sitter dynamics is recovered. In the 
opposite limit (i. e. r 3> +T\), a(r) ~ a x r and the radiation dominance is recovered. In Fig. [5] (plot 
at the left) the exponential damping of the mixing coefficients is numerically illustrated. The curve at 
the top (full line) illustrates the case k — 1. The cases k = 2 and k = 3 are barely distinguishable at 
the bottom of the plot. Notice, always in the right plot, the rather narrow range of times which are 
reported in a linear scale. In the plot at the right the asymptotic values of the mixing coefficients are 
reported for different values of k = k/k max . By fitting the numerical data with with an equation of 
the form given in Eq. ( 12.701) . the value of /3 — 6.33. Different examples can be presented on the same 
line of the one discussed in Fig. [51 While it is clear, from the numerical data, that the decay is indeed 
exponential, the value of /3 may well vary for different models of the transition. The latter observation 
is effectively equivalent to a rescaling of k max for different models of inflation-radiation transition. By 
positing, for instance that /c max — > k mauX /j3 we will have a new k max which differs slightly from k max . It 
is clear that this indetermination on the maximal frequency of the relic graviton spectrum can only 
be solved by endorsing a given model (i.e. by theoretical prejudice) or by having direct measurements 
at those frequencies (which seems to be unlikely in the near future). 
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Figure 6: The spectral energy density of relic gravitons in critical units (plot at the left). The strain 
amplitude is instead reported in the plot at the right. Note that while /iq^gw^, t o) is dimensionless, 
Sh(v,T~o) has dimensions of Hz -1 . The fiducial set of parameters used corresponds to the best fit to 
the WMAP 5-yr data. 

3 Nearly scale- invariant spectra 

The transfer function of the spectral energy density has been numerically computed in the previous 
section and the numerical results have been corroborated by appropriate semi-analytical estimates. 
We are then ready for an explicit calculation of the spectral energy density in the ACDM scenario. In 
subsection 13. li the spectral energy density will be computed in terms of the amplitude transfer function 
and also directly in terms of the transfer function for the spectral energy density. Explicit calculations 
will show that the latter method is more accurate. Subsection 13.21 is devoted to various late time 
effects (e.g. neutrino anisotropic stress, late dominance of dark-energy, progressive diminishment of 
the number of relativistic species) which are certainly present and which affect the amplitude of the 
spectral energy density. Finally, in subsection 13.31 current (and foreseen) sensitivities of wide-band 
interferometers will be briefly compared to the (unfortunately minute) ACDM signal. 



3.1 Spectral energy density in the ACDM paradigm 

In Fig. E] /iq^gw( z/ ; r o) an d the strain amplitude Sh(v,T ) are computed using the transfer function 
for the amplitude discussed (and rederived) in Eqs. (12.281) . (12.291) and (I2.30p . The strain amplitude 
appearing in the plot at the right of Fig. [6] is related to the spectral energy density as in Eq. (12.161) 
which also implies that n GW (j/,r ) can be expressed in terms of Sh(u,To). Indeed, according to Eq. 
(I2.16p . V^k, r) = 4vSh(v, t) and the spectral energy density becomes: 



^gw(^,t) 



4vr 2 



V 3 S h {v,T), 



(3.1) 
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where, in natural units, k = 2irv. By making more explicit the numerical factors and by inverting Eq. 
(13.11) in terms of Shiy, To) we obtain: 



Sh (y, To) 



100 Hz\ 3 
7.981 x l(r ,; ( J /i 2 fi GW (i/,r ) Hz" 1 , 



v 



(3.2) 



where H Q = 3.24078 x 10~ 18 h Hz. The oscillations of Fig. [6] are related to the way the transfer 
function for the amplitude is derived. There are complementary forms of the strategy leading to the 
results of Fig. [6] (see, for instance, |66j). The plots of Fig. [6] have been obtained by using directly 
Eq. (I2.28I) - (I2.30I) inside Eq. (12.431) . This is the procedure used, originally, in [41] (see also [67]). One 
could also define the transfer function as in Eq. (I2.29J) and then, at the level of the spectral energy 
density, compute gk{j) = f' k {r) (see Appendix E of [66J). By working with the transfer function of the 
tensor amplitude the spectral energy density for frequencies v 3> z/ eq is given by: 



h 2 tt GW (v,T ) 



Afh r r — 



v 



fix 



e "max 



7.992 x 10 



-15 



n 



Bn 



0"R0 



di 



0.13267 V4.15 x 10- 5 A 1.4115 x 10 4 Mpc 



(3.3) 
(3.4) 



where g^a^*) is the (comoving) angular diameter distance to decoupling. The dependence upon ^(z*) 
arises because we have to estimate r . In Eqs. ( I3.3p - (l3.4p (as well as in the program used for the 
numerical calculations) there are two complementary options: the first one is to give the angular 
diameter distance to decoupling (which is directly inferred from the CMB data). In the case of the 
5-yr WMAP data alone, d^z*) = 14115 Mpc^ 8 , 8 : this approach has been followed, for instance, in 
|67j . In a complementary perspective it is also possible to take the best fit value of the total matter 
fraction (i.e. ^mo = 0.258 for the case of the WMAP 5-yr data alone) and compute the comoving 
angular diameter distance according to the well know expression for spatially flat Universes: 



d\(z*) = — 



dz 



H ° Jo Jn M0 (i + z) 3 + n A + n R0 (i + z)* 



3.375 



14072 Mpc, 



1090. (3.5) 



The latter strategy has been used, for instance, in [66]. The two strategies are compatible and, 
moreover, this explains why, in Eq. (13. 4p the dependence upon Q M0 does not cancel. In Eq. (13. 3p n T 
denotes, as usual, the tensor spectral index which can be also written as 



n T = -2e+^- In (k/k p ) 



a^ 



r T 



>s " 1) + 



r T 



(3.6) 



If ot = Eq. (I2.23P is recovered and the spectral index is independent on the frequency. In the case 
when ot^O and it is given by Eq. (13.61) the spectral index does depend upon the frequency: in the 
jargon this is often dubbed by saying the the spectral index runs. The frequency-dependent correction 
(i.e. cut) contains the scalar spectral index n s and this is why the value of n s is mentioned in the 
parameters of Fig. The last remark concerning the result of Eq. (13 .4p is that, in the limit, v ^> z/ eq , 
the oscillating terms have been appropriately averaged: this is done by setting the terms going as 
cos 2 (27rz/ro) to 1/2. The latter procedure has been employed, for instance, in the analyses of [671 170] . 
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This procedure is justified in semi- analytical terms but rather odd in a fully numerical context. In 
what follows it will be argued that there is no need of this type of tricks if the spectral energy density 
is computed directly without passing through the transfer function of the tensor amplitude. 

Along this perspective, the results of Fig. [6]should then be compared with Fig. [TJwhere the transfer 
function for the spectral energy density has been consistently employed. In Fig. [7] the spectral energy 
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Figure 7: The spectral energy density of the relic gravitons (plot at the left) and the related Sh(v, r ) 
(plot at the left) for different values of r^ and for the same set of fiducial parameters illustrated in 
Fig. El 

density of the relic gravitons as well as Sh(v, tq) are reported for different values of rx and for the same 
fiducial set of parameters used in Fig. [6j The first salient feature emerging from the comparison of 
Figs. [6] and [7J is that the oscillatory behaviour disappear. The spectra of Fig. [7J have been obtained 
from the direct integration of the mode functions but can be parametrized, according to Eq. fj2.53|) as 



h tt GW (v,r ) 



N n 



AW(z//z/ cq )r T - e 



Vr, 



m 



-2/3- 



4.165 x 10 



-15 



^o^Ro 



(3.7) 



(3.8) 



,4.15 x 10- 5 , 

By comparing Eqs. ( I3.3p -( 13T4"I) to Eqs. (I3.7p -( l3~8j) . the amplitude for v ^> z/ eq differs by a factor which 
is roughly a factor 2. This occurrence is not surprising since Eqs. A3.3J) — fj3.4f) have been obtained 
by averaging over the oscillations (i.e. by replacing cosine squared with 1/2) and by imposing that 
|<7fc| — k\fk\- In Fig. [7] the impact of the variation of tit is also illustrated. Recalling that the WMAP 
5-yr data alone sugggest tt < 0.4, the variation of the spectral energy density is more pronounced than 
the change of the strain power spectrum. This is because of the steepness of Sh(u, r ) in frequency. 



3.2 Anisotropic stress and dark-energy contribution 

The considerations of the previous subsection suggest that the results obtainable with the transfer 
function of the spectral energy density seem to be intrinsically more accurate. The obvious question 
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is of course if we need this precision. There are two answers to this kind of questions. The first 
one is that, of course, the accuracy in the estimate of the ACDM plateau is necessary for comparing 
the theoretical predictions with the data. Therefore it would be strange to treat very accurately the 
tensor contribution to CMB anisotropies but not to wide-band detectors. The second issue is more 
theoretical. In the recent past the community investigated various late time effects which can modify 
the ACDM plateau for v ^> u eq . All these effects compete with the accuracy which is inherent in the 
estimate of the transfer function. This will be the subject of the present subsection. 

Let us therefore start by noticing that, so far, the evolution of the tensor modes has been treated as 
if the anisotropic stress of the fluid was absent. After neutrino decoupling, the neutrinos free stream 
and the effective energy-momentum tensor acquires, to first-order in the amplitude of the plasma 
fluctuations, an anisotropic stress, i.e. 

ST? = -5p8$+Ili flfcn} = II* = 0. (3.9) 

The presence of the anisotropic stress clearly affects the evolution the tensor modes whose evolution 
is then dictated by 

hf + 2Hh( - V 2 h{ = -167iGa 2 U{. (3.10) 

Equation (I3.10p reduces to an integro-differential equation which has been analyzed in [SB] (see also 
[591 ES EI])- The overall effect of collisionless particles is a reduction of the spectral energy density of 
the relic gravitons. Assuming that the only collisionless species in the thermal history of the Universe 
are the neutrinos, the amount of suppression can be parametrized by the function 

F(R U ) = 1 - 0.539^ + 0.134R* (3.11) 

where R v is the fraction of neutrinos in the radiation plasma, i.e. 

r + V " V 3 

In Eq. ( I3.12p N v represents the number of massless neutrino families. In the standard ACDM scenario 
the neutrinos are taken to be massless. In the case R v = (i.e. in the absence of collisionless 
patrticles) there is no suppression. If, on the contrary, R u ^ the suppression can even reach one 
order of magnitude. In the case N v = 3, R v = 0.405 and the suppression of the spectral energy density 
is proportional to JF 2 (0.405) = 0.645. This suppression will be effective for relatively small frequencies 
which are larger than u cq and smaller than the frequency corresponding to the Hubble radius at the 
time of big-bang nucleosynthesis, i.e. z/ bbn of Eq. fll.Sp . 

The effect of neutrino free streaming has been included in Fig. [H] together with the damping effect 
associated with the (present) dominance of the dark energy component. The redshift of A-dominance 
is defined as 

1 + , A ^U(ik) 1/3 . (3,3) 

\a A J Wmo/ 

Consider now the mode which will be denoted as k\, i.e. the mode reentering the Hubble radius at T\. 
By definition k\ = H\a\ must hold. But for r > ta is constant, i.e. H\ = Ho where Hq is the present 
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R v = ——^ r = 0.681 (-^), R^ + R u = l. (3.12) 
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Figure 8: In both plots the contribution of the neutrino free streaming and of the variation in the 
number of degrees of freedom has been taken into account. In the plot at the right the spectral index 
has been allowed to depend upon the frequency. 

value of the Hubble rate. Using now Eq. ( I3.13p . it can be easily shown that k\ = (f2Mo/^A) 1//3 fcH 
where &h = clqHq. The frequency interval between v^ and v^ is rather tiny. Indeed, it turns out that 
^a = ^a/(2vt) is given by 



v K = 2.607 x 10 



-19 



h 







0.719/ V0.258 



n 



MO 



1/3 ( Oa \ 1/3 
0.742 



Hz. 



(3.14) 



For the same choice of parameters of Eq. ( 13.14(1 . vyl = Hq/(2tc) = 3.708 x 10~ 19 Hz which is not 
so different than u\ = 2.607 x 10~ 19 Hz. The adiabatic damping of the mode function across T\ 
reduces the amplitude of the spectral energy density by a factor (fi M0 /Q^) 2 . For the typical choice of 
parameters of Eq. (13.141) we have that the suppression is of the order of 0.12. This class of effects has 
been repeatedly in a number of recent papers [721 E3]- The essence of the effect is captured by the 
following observation. Consider a mode k which reenters before ta- The present value of the amplitude 
Fk{. T ) = fk( T )/ a ( T ) will be adiabatically suppressed since, as repeatedly stressed, in this regime /&(t) 
will simply be plane waves. Consequently, defining as F^ the amplitude at fc* = H*a* when the given 
mode crosses the Hubble radius, we will also have that 



F k (T ) 



OkA /oa\ 



&A / mat V Go / A 



F k =^- 
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k H 



-2 
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M0 



a< 



'it., 



(3.15) 



where the subscripts (in the first equality) denote the time range over which the corresponding redshift 
is computed, i.e. either matter-dominated or A-dominated stages. The second equality follows from 
the first one by appreciating that o(fc*) — rf ~ k~ 2 and by using Eq. (13.131) . Equation (13.151) implies, 
immediately, that the spectral energy density of relic gravitons is corrected in two different fashions. 
For v < z/ H the frequency dependence will be different and will be proportional to f2 GW (z/, r ) oc 
(i / /^'h) TIt ~ 2 (^mo/^a) 2 . Vice versa, in the range v > z/h the frequency dependence will be exactly 
the one already computed but, overall, the amplitude will be smaller by a factor (Qmo/^a) 2 - Two 
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comments are in order. The modification of the frequency dependence is only effective between 10 ! 
0.36 aHz and 0.26 aHz: this effect is therefore unimportant and customarily ignored (see, for instance, 
[63, [72]) for phenomenological purposes. On the other hand, the overall suppression going as (^mo/^a) 2 
must be taken properly into account on the same footing of other sources of suppression of the spectral 
energy density. 

There is, in principle, a third effect which may arise and it has to do with the variation of the 
effective number of relativistic species. The total energy density and the total entropy density of the 
plasma can be written as 

7T 2 27T 2 

Pt = 9 P (T)-T\ s t = g s (T)—T\ (3.16) 

For temperatures much larger than the top quark mass, all the known species of the minimal standard 
model of particle interactions are in local thermal equilibrium, then g p = g s = 106.75. Below, T ~ 175 
GeV the various species start decoupling, the notion of thermal equilibrium is replaced by the notion 
of kinetic equilibrium and the time evolution of the number of relativistic degrees of freedom effectively 
changes the evolution of the Hubble rate. In principle if a given mode k reenters the Hubble radius at 
a temperature T^ the spectral energy density of the relic gravitons is (kinematically) suppressed by a 
factor which can be written as (see, for instance, |72|) 



9p(T k )\ f g s (T k )\ 4 / 3 



9p0 / V 9so 

At the present time g p o = 3.36 and g s o = 3.90. In general terms the effect parametrized by Eq. 
( I3.17P will cause a frequency- dependent suppression, i.e. a further modulation of the spectral energy 
density £Iqw(v, tq). The maximal suppression one can expect can be obtained by inserting into Eq. 
(I3.17P the highest possible number of degrees of freedom. So, in the case of the minimal standard 
model this would imply that the suppression (on f2 GW (z/, r )) will be of the order of 0.38. In popular 
supersymmetric extensions of the minimal standard models g p and g s can be as high as, approximately, 
230. This will bring down the figure given above to 0.29. 

All the effects estimated in the last part of the present section (i.e. free streaming, dark energy, 
evolution of relativistic degrees of freedom) have common features. Both in the case of the neutrinos 
and in the case of the evolution of the relativistic degrees of freedom the potential impact of the effect 
could be larger. For instance, suppose that, in the early Universe, the particle model has many more 
degrees of freedom and many more particles which can free stream, at some epoch. At the same 
time we can say that all the aforementioned effects decrease rather than increasing the spectral energy 
density. Taken singularly, each of the effects will decrease I^gw by l ess than one order of magnitude. 
The net result of the combined effects will then be, roughly, a suppression of f^ow^ r o) which is of 
the order of 3 x 10~ 2 (for 10~ 16 Hz < v < 10~ n Hz) and of the order of 4 x 10~ 2 for u > 10" n Hz. The 
impact of the various damping effects is self-evident by looking at Fig. [HI The effects of the neutrinos 
is visible in the intermediate region where the spectrum exhibits a shallow depression. 



10 We are here enforcing the usual terminology stemming from the powers of 10: aHz (for atto Hz i.e. 10 18 Hz), fHz 
(for femto Hz, i.e. 10~ 15 Hz) and so on. 
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3.3 Sensitivity of wide-band interferometers to ACDM signal 

It is now interesting to compare the ideas discussed in the present section with the sensitivity of 
wide-band interferometers. In Fig. [9] the spectral density of relic gravitons is reported, at a specific 
frequency, as a function of rx- The specific frequency at which Qgw^, r o) is computed is given, as 
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Figure 9: The spectral energy density of the relic gravitons in the context of the ACDM model evaluated 
at the Ligo/Virgo frequency as a function of the tensor-to-scalar ratio. In the plot at the left qt = 
(i.e. the spectral index does not "run" and it is independent upon the frequency. In the plot at the 
right «t is given as in Eq. ( 13. 6ft and the spectral index does depend upon the frequency. 

indicated by v^sr = 100Hz. The subscript LV is a shorthand notation for Ligo/Virgo. If Fig. [9] (plot at 
the left) the tensor spectral index is frequency-independent (i.e. «t = in Eq. (I2.23P ). In the same 
Fig. [9] (plot at the right) «t is allowed to run and ax is given, in terms of the scalar spectral index 
nx as in Eq. (13.61) . 

It is the moment of comparing the theoretical signal with the current sensitivity of wide-band 
interferometers. This figure can be assessed, for instance, from Ref. [22] (see also [3TI I3"6] ) where 
the current limits on the presence of an isotropic background of relic gravitons have been illustrated. 
According to the Ligo collaboration (see Eq. (19) of Ref. [35]) the spectral energy density of a putative 
(isotropic) background of relic gravitons can be parametrized aa 11 !: 



^Gw(^,Tb) = ^ 



GW,/3 



V 



\/3+3 

1 00 Hz ) 



(3.18) 



It is worth mentioning that the parametrization of Eq. ( 13.181) fits very well with Fig. [9] where the pivot 
frequency v L y = 100Hz coincides with the pivot frequency appearing in the parametrization (I3.18p . 
For the scale-invariant case (i.e. (3 = —3 in eq. (13.181) ) the Ligo collaboration sets a 90% upper limit 



n To be completely faithful with the Ligo parametrization the variable (3 will not be changed. It should be borne in 
mind, however, that /3 is used, in the present paper, to quantify the theoretical error on the maximal frequency of the 
relic graviton spectrum (see e.g. Eq. (|2.70p and discussion therein). 
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of 1.20 x 10 -4 on the amplitude appearing in Eq. (I3.18p . i.e. £)gw,-3- Using different sets of data 
(see [34"l [36] ) the Ligo collaboration manages to improve the bound even by a factor 2 getting down 
to 6.5 x 10~ 5 . Keeping an eye on Fig. [9] shows that the current LIGO sensitivity is still too small. 

As far as the ACDM model is concerned, direct detection looks equally hopeless also for the 
advanced interferometers. In the case of an exactly scale invariant spectrum the correlation of the 
two (coaligned) LIGO detectors with central corner stations in Livingston (Lousiana) and in Hanford 
(Washington) might reach a sensitivity to a flat spectrum which is [54"1 1551 [56] 

hi Q GW (v LY , t ) ~ 6.5 x 10- 11 O^P) SNR 2 , u hV = 0.1 kHz (3.19) 

where T denotes the observation time and SNR is the signal to noise ratio. Equation (13.191) is in 
close agreement with the sensitivity of the advanced Ligo apparatus [I] to an exactly scale-invariant 
spectral energy density [921 [931 El] • Equation (13. 19[) together with the plots of Fig. [9] suggest that the 
relic graviton background predicted by the ACDM paradigm is not directly observable by wide-band 
interferometers in their advanced version. The minuteness of /io^Gw( z/ LV; r o) stems directly from the 
assumption that the inflationary phase is suddenly followed by the radiation-dominated phase. 



4 Scaling violations at high frequencies 

According to the results of the previous section, even in the future, the sensitivity of wide-band 
interferometers will be insufficient to reach into the parameter space of the ACDM scenario. The 
accurate techniques introduced in the present paper seem therefore a bit pleonastic. In this section 
the opposite will be argued insofar as the spectral energy density of relic gravitons may well be 
increasing (rather than decreasing) as a function of the frequency v. 

The late and early time effects conspire, in the ACDM paradigm to make the spectral energy 
density slightly decreasing at high frequencies (see e.g. Fig. [S]). Different thermal histories allow for 
scaling violations which may also go in the opposite direction and make the spectral energy density 
increasing (rather than decreasing as in the ACDM case) for typical frequencies larger than a pivotal 
frequency z/ s which is related to the total duration of the stiff phase. If the stiff phase takes place 
before BBN, then z/ s > 10 -2 nHz. If the stiff phase takes place for equivalent temperatures larger than 
100 GeV, then v s > //Hz. Finally, if the stiff phase takes place for T > 100 TeV, then z/ s > mHz. 

In the early Universe, the dominant energy condition might be violated and this observation will 
also produce scaling violations in the spectral energy density [71]. If we assume the validity of the 
ACDM paradigm, a violation of the dominant energy condition implies that, during an early stage of 
the life of the Universe, the effective enthalpy density of the sources driving the geometry was negative 
and this may happen in the presence of bulk viscous stresses [71] (see also [751 [75] for interesting 
reprises of this idea). In what follows the focus will be on the more mundane possibility that the 
thermal history of the plasma includes a phase where the speed of sound was close to the speed of 
light. In Eqs. (12.541) and (12.551) the stiff evolution has been parametrized in terms of an effective (i.e. 
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fluid) description which can be realized in diverse models not necessarily related to a fluid behaviour. 
If the energy-momentum tensor of the sources of the geometry is provided by a scalar degree of freedom 
(be it for instance ip) The effective energy density, pressure and anisotropic stress of tp will then be, 
respectively, 

p* = ( € + v) + i(^) 2 , Pv = (t.-v\- i-Jd^) 2 , (4.i) 



n?(v) = -- 



2 ) 2a 2V ™' ' * v V2 ) 6a 2 

(4.2) 



d i (pd j (p--{d k (p) 2 5 J i 



Equation (14. ip imply that the effective barotropic index for the scalar system under discussion is simply 
given by 



Pv 



P * (% + v)+&(d k <p)* 



(4.3) 



If <p 2 3> V and tp 2 3> (dk(p) 2 /a 2 , then p v ~ p v : in this regime the scalar field behaves as a stiff fluid. If 
F>^) 2 > (dk<f) 2 /a 2 , then w v ~ — 1: in this regime the scalar field is an inflaton candidate. Finally 
if (dk<p) 2 /a 2 ^> tp 2 and (dk<p) 2 /a 2 ^> V, then w v ~ —1/3: in this regime the system is gradient- 
dominated. Of course also intermediate situations are possible (or plausible). 

From the purely phenomenological point of view it is not forbidden (by any phenomenological 
consideration) to have a sufficiently long stiff phase. This was the point of view invoked in [ST] (see 
also [521 [S3]) where it was also suggested that the spectral energy density of relic gravitons may increase 
with frequency. The presence of a phase dominated by the kinetic energy of a scalar degree of freedom 
(typical of quintessence models) became more compelling also in the light of the formulation of the 
so-called quintessential inflationary models jSH] where the inflaton field practically does not decay and 
it is identified with the quintessence field. If there is some delay between the end of inflation and the 
onset of radiation the maximal wavenumber of the spectrum will be given by: 

where a = 2/[3(w t + 1)] is related to the specific kind of stiff dynamics (indeed, w t > 1/3). Equation 
(14.41) can also be written as 

A; max = MpE-^lp) 172 ^*), v M = WW- (4.5) 

where 

/ H \ a ~ 1 / H r \ 1 / 2_Q 

In the case S = 0(1) (as it happens in the case a = 1/3 if the initial radiation is in the form of 
quantum fluctuations) vm — 100 GHz, more precisely: 

z/ max = 1.177 x lO 11 ^ 1 ( k2 ° nR0 - ) V4 Hz. (4.7) 

V4.15 x 10-V 
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The strategy will now be to parametrize the violation of scale invariance in terms of the least possible 
number of parameters, i.e the frequency u B (defining the region of the spectrum at which the scaling 
violations take place) and the slope of the spectrum arising during the stiff phase. Of course the 
frequency u s can be dynamically related to the frequency of the maximum and, consequently, the first 
parameter can be trated for S. The slope of the spectrum during the stiff phase depends upon the 
total barotropic index and can therefore be traded for w t . Assuming the presence of a single stiff 
(post-inflationary) phase we will have that 

where the second equality follows from the first by using the relation of H T to X dictated by Eq. (j4.8p . 
From Eq. (14. 8 j) the frequency turns out to be: 

u a = 1.173 x lO^/d-sc*) {7TeAn) w^ ( 4i ^ R i ° Q _ 5 ) 1/4 Hz. (4.9) 

The quantity S is always smaller than 1 or, at most, of order 1. This is what happens within specific 
models. For instance, if the radiation present at the end of inflation comes from amplified quantum 
fluctuations (i.e. Gibbons-Hawking radiation), quite generically, at the end of inflation p r ~ H 4 . More 
specifically 

30" *""" 48071-* 

In Eq. (14.101) N e g is the number of species contributing to the quantum fluctuations during the quasi- 
de Sitter stage of expansion. In [61] (see also [521 [53], [59]) it has been argued that this quantity could 
be evaluated using a perturbative expansion valid in the limit of quasi-conformal coupling. It should 
be clear that N eT j is conceptually different from the number of relativistic degrees of freedom g p . Given 
H and N e s the length of the stiff phase is fixed, in this case, by [59J 



7T 4 N eS H 

Pr = ^N eS T A H = -^—. (4.10) 



<i\ 



i 



a .x 3(™+i) , 2 faA 2 / a 



A# 4 P =H 2 MU^) =H 2 Ml[^\ , (4.11) 

V a T J \a T J \a r J 

where we used the fact that a = 2/[3(w + 1)] and where we defined A = iV c ff/(4807r 2 ). Equation ( 14.111) 
implies that 

Oi \ , a / H \ T^ ( H T \ , 1 / H \ (1-2q) 

Using the second relation in Eq. (14. 12j) and Eq. (14. 6p . it turns out that S = A 1//4 , which is always 
smaller than 1 and, at most, 0(1). Instead of endorsing an explicit model by pretending to know the 
whole thermal history of the Universe in reasonable detail, it is more productive to keep S as a free 
parameter and to require that the scaling violations in the spectral energy density will take place before 
BBN. The variation of S, w and r^ can be simultaneously bounded. The essential constraint which 
must be enforced in any model of scaling violations implies that the frequency v s must necessarily 
exceed z/ bbn (see Eq. (II. 8p ). This requirement guarantees that the stiff dynamics will be over by the 
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Figure 10: In the plot at the left the bounds on u B are illustrated in terms of w t (i.e. the barotropic 
index of the stiff phase). The quantity £ (defined in Eq. (14.61) ) depends, both, on the Hubble rate at 
the end of inflation and on the Hubble rate at the onset of radiation. At the right the exclusion region 
is phrased in terms of rx (i.e. the tensor to scalar ratio) for different values of w t . The shaded areas 
represent the excluded regions. 

time light nuclei start being formed. In a complementary approach one might also require that u s > u cw 
where v ew corresponds to the value of the Hubble rate at the electroweak epoch, i.e. 



3.998 x 10" 



9p 



i/i 



T 



hlQ. 



0"R0 



1/1 



Hz. 



.106.757 V100 GeVy V4.15 x 10- 5 , 
Finally, yet a different requirement could be to impose that v > z/t cv where z/Tev is defined as 



z/ TeV = 4.819 x 10 



-3 



9 P 



1/4 



T, 



tin 



0"R0 



1/1 



Hz. 



(4.13) 



(4.14) 



.228.757 V100 TeV7 V4.15 x 10- 5 . 

The condition v > u^ev (as opposed to v > u cw ) would imply that the stiff age did already finish by 
the time the Universe had a temperature of the order of 100 TeV when, presumably, the number of 
relativistic degrees of freedom was much larger than in the minimal standard model l 12 l . 

The constraints on E, w t and r^ are summarized in Fig. [10] The value of E controls the position of 
the frequency at which the nearly scale-invariant slope of the spectrum will be violated. The barotropic 
index w t is taken to be always larger than 1/3 (by definition of stiff fluid) and with a maximal value 
of 1. In Fig. [TD] values of w t as large as 2 have been allowed just for completeness since some authors 
like to speculate that models with w t > 1 do not violate causality constraints. It is amusing to notice 
that the cases w t > 1 are, anyway, totally irrelevant from the phenomenological point of view. In these 
cases, in fact, the detectability prospects are forlorn (see section [5]). 

In Fig. [10] (plot at the right) the different curves denote, respectively, the cases z/ s = u hhn (full 
line), v s = v cw (dashed line) and z/ s = ^TeV (dot-dashed line). To be compatible with the corresponding 



12 In Eq. (14. 14[) the typical value of g p is the one arising in the minimal supersymmetric extension of the standard 
model. 
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Figure 11: The spectral energy density of the relic gravitons coming from the stiff ages. In the plot at 
the left Wt — 1 while in the plot at the right w t = 0.6. In both plots the value of E has been fixed to 
0.15. 

constraint we have to be above each curve and the shaded region are excluded. Of course since different 
curves are present, we decided to shade the region which might correspond, according to theoretical 
prejudice to the most typical choice of parameters. 



5 Relic gravitons and from the stiff age 

5.1 Spectral energy density in the minimal TACDM scenario 

The conclusion of the previous section has been that it is indeed plausible to parametrize the scaling 
violations (at high frequency) in terms of two parameters, i.e. the typical frequency u s at which 
scaling violations occur and the typical slope of the spectral energy density for v > u 8 . The latter 
framework has been dubbed TACDM for tensor-ACDM [77J. The two supplementary parameters 
physically depend upon the sound speed during the stiff phase (i.e. c st ) and the threshold frequency 
(i.e. z/ s ). Besides c st and u s , there will also be r^ which controls, at once, the normalization and the 
slope of the low-frequency branch of the spectral energy density. The remaining six parameters of the 
underlying ACDM model will be fixed, just for illustration, to the best fit of the 5-yr WMAP data 
alone [3, El [TJ El [9] . In the numerical program used to compute the spectral energy density of the relic 
gravitons the putative values of the cosmological parameters can be changed at wish. 

In spite of their present sensitivities [35j (see also Eq. (13.181) and discussion therein) terrestrial 
interferometers might be able, one day, to provide a prima facie evidence of relic gravitons. The present 
numerical approach will then be instrumental not only in setting upper limits but also in providing 
global fits of the cosmological observables within the TACDM model. According to this perspective, in 
the future the three (now available) cosmological data sets will be complemented by the observations 
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of the wide-band interferometers. Therefore, different choices of cosmological parameters (like, for 
instance, the various critical fractions of matter and dark energy) could slightly change the typical 
frequencies of the relic graviton spectrum as well as other features in the low-frequency region of the 
spectral energy density. 

In the absence of any tensor contribution (i.e. rx = 0) the 5-yr WMAP data alone imply: 

(fibo, tt c0 , tt A , h , n s , t) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087). (5.1) 

If the tensors are included (i.e. tt 7^ 0) but without any running of the scalar spectral index the 
parameters (inferred from the 5-yr best fit to the WMAP data alone) slightly change and become: 

(ftbo, n c0 , Sl A , ho, n s , r,r T ) = (0.0417, 0.188, 0.770, 0.751, 0.986, 0.090, < 0.43). (5.2) 

In Eq. (15.2J1 the last entry of the array contains rx and it is actually an upper limit (95% CL) 
corresponding to the first row appearing in Tab. 11.21 Various other examples could be provided by 
considering, for instance, the combinations listed in Tab. 11.21 In the numerical examples reported 
here, the ACDM parameters will be fixed to their best fit values as they are reported in Eq. (15.11) . In 
this situation the tensor contribution will be parametrized not only by r-r, but also by c s t and v s . The 
bounds on r T are spelled out in Tab. 11.21 

In both plots of Fig. [TTJthe parameters are fixed to the values reported in Eq. f)5.ip . In Fig. [TT] 
(plot at the left) the ACDM scenario is complemented by a stiff phase with w t = 1 and for different 
values of r-r- Always in Fig. [TTJthe value of the barotropic index is slightly reduced from ltow t = 0.6. 
In both plots of Fig. [Til «t/0 and its valudjfl is given by Eq. (13.61) . The effect associated with 
a slight frequency variation of the tensor spectral index is rather modest so that it can be hardly 
distinguished from «t = except when rx is sufficiently large. A similar occurrence can be observed 
in the two plots reported in Fig. [9] 

The infrared branch of the spectrum in both plots of Fig. [TT] reproduces the results of Fig. [HI 

As soon as the frequency increases from the aHz up to the nHz (and even larger) the spectral 
energy density increases sharply in comparison with the nearly scale-invariant case (see, e.g. Figs. [8] 
and [9]) where the spectral energy density was, for v > nHz, at most 0(1O~ 16 ). In the case of Fig. [TTJ 
the spectral energy density is clearly much larger. The accuracy in the determination of the infra-red 
branch of the spectrum is a condition for the correctness of the estimate of the spectral energy density 
of the high-frequency branch. The plots of Fig. [TT] demonstrate that the low-frequency bounds on r^ 
do not forbid a larger signal at higher frequencies. 

A decrease of r^ implies a suppression of the nearly scale-invariant plateau in the region u eq < 
v < v B . At the same time the amplitude of the spectral energy density still increases for frequencies 
larger than the frequency of the elbow (i.e. u s ). The latter trend can be simply understood since, 
at high frequency, the transfer function for the spectral energy density grows faster than the power 



13 If not otherwise specified, the value of the scalar spectral index used to compute ot is consistent with the 5-yr best 
fit to the WMAP data alone. 
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spectrum of inflationary origin. For instance, in the case w t = 1 and neglecting logarithmic corrections, 
Ocw( i/ ) r o) ot z/™ T+1 for v ^> z/ s . Now, recall that n^ is given by Eq. (13.61) . If r T — > 0, the combination 
(«x + 1) will be much closer to 1 than in the case when, say, r^ ~ 0.3. This aspect can be observed 
in both plots of Fig. (TT] where different values of tt have been reported. By decreasing the u>t from 1 
to, say, 0.6 the extension of the nearly flat plateau gets narrower. This is also a general effect which is 
particularly evident by comparing the two plots of Fig. [Til The slope of the high-frequency branch of 
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Figure 12: The bounds stemming from the amount of extra-relativistic species at the epoch of the 
synthesis of light nuclei are applied to the relic graviton spectra from the stiff epoch. As indicated 
the parameters of the underlying ACDM model are fixed to the best fit derived from the 5-yr WMAP 
data alone (see also Eq. (15.11) ). 

the graviton energy spectrum can be easily deduced with analytic methods and it turns out to to be 

d In Ogw 6u> t — 2 



dlnu 



3w t + 1 ' 



v > i/ s , 



(5.3) 



up to logarithmic corrections. The result of Eq. (15.31) stems from the simultaneous integration of the 
background evolution equations and of the tensor mode functions according to the techniques described 
in section[3J The semi-analytic estimate of the slope (see [SI]) agrees with the results obtained by means 
of the transfer function of the spectral energy density. In Fig. 0] (plot at the left), for k = k/k s > 1 
T?(k) — K which is consistent with Eq. ( I5.3P in the case w t = 1. The logarithmic corrections arising 
in the case w t = 1 (see, for instance, Eq. (12.571) ) have a simple analytic interpretation which is evident 
from the results reported in Eqs. (I2.65p - (|2.66p and (I2.67I) - (I2.68I) for the mixing coefficients in the case 
w t = 1. 

5.2 Phenomenological constraints 

The spectra illustrated in Fig. [TTJ (as all the spectra stemming from the stiff ages) must be compatible 
not only with the CMB constraints (bounding, from above, the value of rx) but also with other two 
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classes of constraints, i.e. the pulsar timing constraints [78j [79] and the big-bang nucleosynthesis 
constraints [80j EH [82] . The pulsar timing constraint demands 

fi(^ P uisar,i- ) < 1.9 x 10~ 8 , z/ pu i sar ~ 10 nHz, (5.4) 

where pulsar roughly corresponds to the inverse of the observation time along which the pulsars timing 
has been monitored. Assuming the maximal growth of the spectral energy density and the minimal 
value of u s , i.e. z/ bbn we will have 

hltt GW (v, T ) OC U, V>V S ^ ffatan- (5-5) 

Since f pu isar — lO^bbn, Eq. (15.51) implies that /iq^qw (^pulsar, To) — 10 -13 or even 10~ 14 depending 
upon rx- But this value is always much smaller than the constraint stemming from pulsar timing 
measurements. If either u s ^> z/ bbn or c st < 1 the value of ^o^gw (^pulsar, To) will be even smalleii 14 !. 
Consequently, even in the extreme cases when the frequency of the elbow is close to z/ bbn , the spectral 
energy density is always much smaller than the requirement of Eq. (15. 41) . The conclusion is that the 
pulsar timing bound is not constraining for the TACDM model. 



It is well known that the most significant constraint on the stiff spectra stems from BBN [49, 
W\\ . Being massless, gravitons can increase the expansion rate at the BBN epoch. To avoid the 
overproduction of 4 He, the number of relativistic species must be bounded from above. The BBN 
bound is customarily expressed in terms of (equivalent) extra fermionic species. According to Eq. 
(I3.16p . during the radiation-dominated era, the energy density of the plasma can be written as p t = 
g p (7r 2 /30)T 4 where T denotes here the common (thermodynamic) temperature. An (ultra)relativistic 
fermion species with two internal degrees of freedom and in thermal equilibrium contributes 2 • 7/8 = 
7/4 = 1.75 to g p . Before neutrino decoupling the contributing relativistic particles are photons, 
electrons, positrons, and N u = 3 species of neutrinos, giving g p = 10.75. 

The neutrinos have decoupled before electron-positron annihilation so that they do not contribute 
to the entropy released in the annihilation. While they are relativistic, the neutrinos still retain 
an equilibrium energy distribution, but after the annihilation their (kinetic) temperature is lower, 
T v = (4/ll) 1 ' 3 T. Thus g p = 3.36 after electron-positron annihilation. By now assuming that there 
are some additional relativistic degrees of freedom, which also have decoupled by the time of electron- 
positron annihilation, or just some additional component p x to the energy density with a radiation-like 
equation of state (i.e. p x = p x /3), the effect on the expansion rate will be the same as that of having 
some (perhaps a fractional number of) additional neutrino species. Thus its contribution can be 
represented by replacing N u with N v + AA;, in the above. Before electron-positron annihilation we 
have p x = (7/8)AA„p 7 and after electron-positron annihilation we have p x = (7/8)(4/ll) 4 / 3 AA„p 7 ~ 
0.227 AA^p 7 . 

The critical fraction of CMB photons can be directly computed from the value of the CMB tem- 
perature and it is notoriously given by tiffin = p 7 /p cr i t = 2.47 x 10~ 5 . If the extra energy density 



14 This conclusion follows immediately from the hierarchy between ^ pu i sar and t'bbn- If either c st < 1 or i/ s > ^bbm 
^o^gw can only grow very little and certainly much less than required to violate the bound of Eq. ([57 
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hln x = tf^ = 5.61 x 10-*AN U ( - ££%-r ) • (5.6) 



hl£™ X n GW (v,T )d\nis = 5M x 10 " 6A ^( 2 4 ^ -5 )> ( 5 - 7 ) 



component has stayed radiation-like until today, its ratio to the critical density, Q x , is given by 

p c ""'" ' '" ""\2A7x 10-6 
If the additional species are relic gravitons, then [EDI EH E2] : 

rv *>°* ... „ / /ln^- 

'^bbn 

where z/ bbn and z/ max are given, respectively, by Eqs. (11.81) and (14.71) . Thus the constraint of Eq. (15.71) 
arises from the simple consideration that new massless particles could eventually increase the expansion 
rate at the epoch of BBN. The extra-relativistic species do not have to be, however, fermionic [81] 
and therefore the bounds on AN V can be translated into bounds on the energy density of the relic 
gravitons. 

A review of the constraints on AN V can be found in [81] . Depending on the combined data sets 
(i.e. various light elements abundances and different combinations of CMB observations), the standard 
BBN scenario implies that the bounds on AN U range from AN V < 0.2 to AN U < 1. Similar figures, 
depending on the priors of the analysis, have been obtained in a more recent analysis [82]. All the 
relativistic species present inside the Hubble radius at the BBN contribute to the potential increase 
in the expansion rate and this explains why the integral in Eq. (15.71) must be performed from z^bbn to 
^max (see also [55] where this point was stressed in the framework of a specific model). 

The existence of the exponential suppression for v > z/ max (see Fig. [IT]) guarantees the convergence 
of the integral also in the case when the integration is performed up to v — > oo. The constraint of 
Eq. (15 .7p can be relaxed in some non-standard nucleosynthesis scenarios [81] , but, in what follows, the 
validity of Eq. (I5.7P will be enforced by adopting AN U ~ 1 which implies, effectively 

hlJ^tt GW (v,T )dlnv < 5.61 x 10- 6 ( 24 ^ _ 5 ). (5.8) 

The models illustrated in Fig. [TTJ are on the verge of saturating the bounds of Eqs. (J5.7p - (j5.8p . 
This conclusion stems directly from the form of spectral energy density: the broad spike dominates 
the (total) energy density of relic gravitons which are inside the Hubble radius at the time of big bang 
nucleosynthesis. A practical way of enforcing the bounds of Eqs. (15.7]) and (15. 8p is to integrate around 
the maximum of the curves depicted in Fig. [TTJ In Fig. [T2] the energy density of the relic gravitons 
inside the Hubble radius at the nucleosynthesis epoch is reported in the case t"t = 0.1 and for different 
values of E. In the plot at the left n s = 0.963 as implied by the WMAP 5-yr data alone. In the plot 
at the right n s = 0.8. The two horizontal lines illustrate the bounds of Eqs. (J5.7p ~ (j5.8p in the cases 
AN U = 1 (full line) and AN U = 0.2 (dashed line) which are, respectively, the least constraining and 
the most constraining situations contemplated by current analyses. In both cases the allowed region 
of the parameter space stays below the horizontal lines. As the scalar spectral index diminishes, the 
constraints are better satisfied since n s controls «t and, consequently, the frequency dependence of 
the tensor spectral index tit (see Eq. ( 13. 6p ) in the case «t ^ 0. 
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Figure 13: The detectability constraints (full lines in both plots) stemming from the putative sensitiv- 
ities of wide-band interferometers in their advanced version. The points corresponding to the spectral 
energy density should lie above the full lines to be potentially interesting for those instruments. 
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Figure 14: The spectral energy density is illustrated for small values of w t and different values of S 
(plot at the left). In the plot at the right S = 0.2 and w t = 0.6. 

5.3 Detectability prospects 

The results presented in the previous subsection suggest that if r^ is bounded from above by the 
cosmological data sets (see e.g. Tab. 11.21) . a detectable signal is expected for vysi — 100 Hz for 
0.35 < w t < 0.61. In this case, following the parametrization of the Ligo collaboration we could say 
that the expected signal can be parametrized as 



^Gw( z/ ? r o) — ^GW 



V 



100 Hz 



n T + 



6mt- 2 
3m t + l 



(5.9) 



which mirrors Eq. (19) of Ref. [35] where the upper limits on the amplitude Qqw have been set. This 
range turns out to be compatible with the bounds of Eqs. (I5.7p - (I5.8I) . The relation (15. 9p could be 
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used by the experimenters to set bounds on Qqw m the same way as upper bounds are obtained in 
the case of nearly scale-invariant spectra (see Eq. (I3.18P and discussion therein). Clearly, from the 
theoretical point of view, flew changes by varying the various TACDM parameters. 

For instance, by lowering w t , /iq^gw(^) r o) increases for v = z/ L v — 0.1 kHz. This trend can 
be inferred from Fig. [13] where the spectral energy density is evaluated exactly for v = v-^y. To 
be detectable by wide band interferometers the parameters of the TACDM must lie above the full 
lines. The region of low barotropic indices emerging neatly from Fig. [13], leads to spectral energy 
densities which are progressively flattening as w t diminishes towards 1/3. Low values of w t bring the 
frequency of the elbow, i.e. u s below 10~ 10 Hz which is unacceptable since it would mean that, during 
nucleosynthesis, the Universe was dominated by the stiff fluid. In Fig. [10] (plot at the left) the region 
above the full line corresponds to a range of parameters for which v s > u^ n : in such a range a decrease 
of w t demands an increase of S. 

The occurrence described in the previous paragraph is illustrated in Fig. [141 where, at the left, w t = 
0.5 and the values of E are the same ones illustrated in Fig. [131 The full, dashed and dot-dashed curves 
illustrated in Fig. [T3] (plot at the left) are incompatible with phenomenological considerations since 
the frequency of the elbow is systematically smaller than z/btm- Once more, this choice of parameters 
would contradict the bounds of Fig. CG)] and would imply that the stiff ] phase is not yet finished at 
the BBN time. In the left plot of Fig. [TH the diamonds denote a model which is compatible with 
BBN considerations but whose signal at the frequency of interferometers is rather small (always three 
orders of magnitude larger than in the case of conventional inflationary models). 

The compatibility with the phenomenological constraints demands that the parameters of the 
TACDM paradigm must lie above the full lines of Fig. [10] The requirements of Fig. [10] suggest, 
therefore, that £ should be raised a bit. In this case the frequency of the elbow gets shifted to the 
right but, at the same time, the overall amplitude of the spike diminishes. The putative amplitude 
remains still much larger than the conventional inflationary signal reported in Fig. [9] 

In Fig. [14] (plot at the right) £ = 0.2 and w t = 0.6. The tensor spectral index is allowed to depend 
upon frequency according to Eq. (13.61) (i.e. «t ^ 0). Two different values of n s are reported. In 
the example of Fig. [TH the phenomenological bounds are all satisfied. In Fig. [15] the spectral energy 
density of the relic gravitons is illustrated as a function of rx for a choice of parameters which is 
compatible with all the bounds applicable to the stochastic backgrounds of the relic gravitons. The 
three curves refer to three different frequencies, i.e. 0.1 kHz, 1 kHz and 10 kHz. Indeed, if the spectrum 
is nearly scale-invariant (as in the case o Fig. [9]) we can compare the potential signal with the central 
frequency of the window. If the signal increases with frequency it is interesting to plot the same curve 
for some significant frequencies inside the window of wide-band interferometers. Even if the frequency 
window extends from few Hz to 10 kHz the maximal sensitivity is in the central region and depends 
upon various important factors which will now be briefly discussed. 

To illustrate more quantitatively this point we remind the expression of the signal-to-noise ratio 
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Figure 15: The graviton energy spectrum is illustrated, in the TACDM scenario, for v = vysi and as a 
function of rx- As in Fig. [9] at the left ot = while, at the right, ax 7^ 0. 



(SNR) in the context of optimal processing required for the detection of stochastic backgrounds: 

7 V)^ W (Z,,T ) 1 1/2 



SNR 2 



ml 




du 



.(1) 



^S^(u)Sk 



(2), 



(5.10) 



(F depends upon the geometry of the two detectors and in the case of the correlation between two 
interferometers F = 2/5; T is the observation time). In Eq. (15.1 Op . S^ k \f) is the (one-sided) noise 
power spectrum (NPS) of the fc-th (k = 1, 2) detector. The NPS contains the important informations 
concerning the noise sources (in broad terms seismic, thermal and shot noises) while 7(1/) is the 
overlap reduction function which is determined by the relative locations and orientations of the two 
detectors. In [51] Eq. (15. 10)) has been used to assess the detectability prospects of gravitons coming 
from a specific model of stiff evolution with w t = 1. At that time the various suppressions of the low- 
frequency amplitude as well as the free-streaming effects were not taken into account. Furthermore, 
the evaluation of the energy transfer function was obtained, in [56], not numerically but by matching 
of the relevant solutions. We do know, by direct comparison, that such a procedure is justified but 
intrinsically less accurate than the one proposed here. It would be interesting to apply Eq. (15.101) 
for the (more accurate) assessment of the sensitivities of different instruments to a potential signal 
stemming from the stiff age l 15 l . 



Equation (15.101) assumes that the intrinsic noises of the detectors are stationary, Gaussian, un- 
correlated, much larger in amplitude than the gravitational strain, and statistically independent on 



15 For intermediate frequencies the integral of Eq. (|5.10|) is sensitive to the form of the overlap reduction function 
which depends upon the mutual position and relative orientations of the interferometers. The function 7(1/) effectively 
cuts-off the integral which defines the signal to noise ratio for a typical frequency v ~ l/(2d) where d is the separation 
between the two detectors. Since £Iqw increases with frequency (at least in the case of relic gravitons from stiff ages) 
at most as v and since there is a v~ & in the denominator, the main contribution to the integral should occur for v < 0.1 
kHz. This argument can be explicit verified in the case of the calculations carried on in [51] and it would be interesting 
to check it also in our improved framework. 
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the strain itself [92j (93J [93] . The integral appearing in Eq. (j5.10p extends over all the frequencies. 
However, the noise power spectra of the detectors are defined in a frequency interval ranging from 
few Hz to 10 kHz. In the latter window, for very small frequencies the seismic disturbances are the 
dominant source of noise. For intermediate and high frequencies the dominant sources of noise are, 
respectively, thermal and electronic (i.e. shot) noises. The wideness of the band is very important 
when cross-correlating two detectors: typically the minimal detectable /iq^gw wm become smaller 
(i.e. the sensitivity will increase) by a factor l/yAuT where Av is the bandwidth and T, as already 
mentioned, is the observation time. Naively, if the minimal detectable signal (by one detector ) is 
/iqO gw ~ 10~ 5 , then the cross-correlation of two identical detector with overlap reduction j(y) = 1 
will detect /^gw ^ 10~ 10 provided Av ~ 100 Hz and T ~ O(lyr) (recall that lyr = 3.15 x 10 7 Hz _1 ). 
The achievable sensitivity of a pair of wide band interferometers crucially depends upon the spectral 
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Figure 16: The graviton energy spectrum is illustrated, in the TACDM scenario, for v = vq and as a 
function of r T . As in Figs. 191 and [151 at the left «t = while, at the right, «t ^ 0. 



slope of the theoretical energy spectrum in the operating window of the detectors. So, a flat spectrum 
will lead to an experimental sensitivity which might not be similar to the sensitivity achievable in 
the case of a blue or violet spectra. Previous calculations [H31 [55J [56] showed that, however, to get 
a reasonable idea of the potential signal it is sufficient to compare the signal with the sensitivity to 
flat spectrum which has been reported in Eq. (I3.19p . Of course any experimental improvement in 
comparison with the values of Eq. ( 13. 19ft will widen the detectability region by making the prospects 
of the whole discussion more rosy. 

In the TACDM paradigm the maximal signal occurs in a frequency region between the MHz and 
the GHz. This intriguing aspect led to the suggestion [531 EH] that microwave cavities [83J can be 
used as GW detectors precisely in the mentioned frequency range. Prototypes of these detectors [S3] 
have been described and the possibility of further improvements in their sensitivity received recently 
attention [851 EH [87J [HU [891 EQ]. Different groups are now concerned with high-frequency gravitons. 
In [86] the ideas put forward in [831 EH EH] have been developed by using electromagnetic cavities 
(i.e. static electromagnetic fields). In [HTJ, ESI ES] dynamical electromagnetic fields (i.e. wave guides) 



41 



have been studied always for the purpose of detecting relic gravitons. In [89] an interesting prototype 
detector was described with frequency of operation of the order of 100 MHz (see also [91J). In Fig. 
[TBI the value of the spectral energy density is reported for v = uc where vq defines the frequency of 
operation of a given electromagnetic detector. In Fig. [16] vq is taken in the MHz range. In both plots 
the horizontal lines denote the bounds of Eqs. ( 15.71) and ( 15. 8ft for two typical values of AN U (i.e., more 
specifically, AN U = 1 and AN U = 0.2). To be compatible with the bounds the values of the spectral 
energy density must be smaller than the horizontal lines. The region of large r^ (i.e. r^ ~ 0{1)) 
is already excluded from CMB upper limits: the plots have been extended also in that region for 
comparison with the analog plots illustrated in Fig. [9] 

Absent direct tests on the thermal history of the plasma prior to neutrino decoupling, the current 
bounds on a tensor component affecting the initial conditions of the CMB anisotropies (and polariza- 
tion) do not forbid a potentially detectable signal for typical frequencies compatible with the window 
of wide-band interferometers. The numerical approach described in the present paper allows for a 
sufficiently accurate estimate of the spectral energy density of the relic gravitons. In the context of the 
class of models analyzed here it is plausible to imagine, in the years to come, a rather intriguing synergy 
between large-scale observations (e.g. CMB physics, measurements of the matter power spectrum and 
supernovae) and small scale observations such as the ones conducted by wide-band interferometers in 
the range between few Hz and 10 kHz. 
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